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Absolute  emission  coefficient  measurements  on  arc-generated 

uranium  plasmas  in  local  thermodynamic  equilibrium  are  described  for 

0 

a  wavelength  bandwidth  of  1050  to  6000  A.  Low-  and  high-pressure  arcs 
were  investigated  for  their  emission  properties,  characteristic  tem- 
peratures and  uranium  partial  pressures.  Temperatures  from  5500  to 
8000  K  and  uranium  partial  pressures  from  0.001  to  0.01  atm  were  found 
at  the  arc  centerline.  The  new  emission  data  are  compared  with  other 
similar  experimental  results  and  to  existing  theoretical  calculations. 
The  effects  of  cold-layer  UF^  photoabsorption  on  uranium  plasma  emis- 
sion characteristics  are  established  for  UF^  molecular  densities 

6 
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ranging  from  1.0  x   10       tol.OxlO       cm"     and   layer  thickness   from 
1.0  to  5.0  cm. 


VI 


Ah   iH-itio   atomic  structure  calculations  were  made  using  relativistic 
Hartree  exchange  wavefunctions ,  from  which  oscillator  strength  distribu- 
tions were  computed  for  transition  arrays  of  interest.  These  calcula- 
tions give  supporting  evidence  as  to  the  credibility  of  the  measured 
emission  at  various  wavelengths,  particularly  in  the  vacuum  ultraviolet. 
It  is  suggested  that  a  consistent  picture  as  to  the  nature  of  uranium 
plasma  emission,  at  these  plasma  conditions,  emerges  and  the  capability 
now  exists  to  successfully  compute  major  emission  features  of  uranium 
and  other  complex  atomic  systems. 
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I,        INTRODUCTION 

With  the  advent  of  the  nuclear  age  and  the  subsequent  strong  trend 
toward  development  of  uranium-based  technology,  research  priorities  con- 
cerning the  nature  of  the  uranium  atom  have  acquired  substantially  in- 
creased importance.  Initial  involvement  in  metallurgical  and  nucleonic 
properties  vjas  largely  due  to  the  apparent  need  for  weapon  fabrication 
and  later  some  aspects  of  reactor  technology.  By  1947  a  new  area  of 
interest  was  the  study  of  uranium  plasmas  at  high  temperature  (6  keV), 
for  a  physical  understanding  of  energy  release  caused  by  nuclear  detona- 
tions. This  was  the  first  significant  attempt  to  model  the  uranium  atom 
using  a. quantum  and  statistical  mechanical  basis  from  which  was  extracted 
thermodynamic  and  optical  information.  Eventually,  interest  was 
generated  in  lower  temperature  (0.5  -^  5  eV)  uranium  plasmas  because  of 

the  potential  usefulness  of  plasma  core  reactors  as  a  means  for  space 

2-4 
propulsion  and  possibly  as  an  energy  source  for  MHD  power  generation. 

More  recently,  v/orldwide  need  for  uranium  isotope  enrichment  using  laser 

5  1  2 
processes    has  opened  a  new  and  significant  area  which  is  stimulating 

much  basic  research  of  the  uranium  atom.  Uranium  plasma  research  is  also 

influencing  the  development  of  nuclear-pumped  laser  systems.  This 

■]  o   "1-1 

particular  program  has  had  some  recent  breakthroughs  '   which  will 
likely  increase  the  research  momentum  on  the  study  of  neutral  and  once- 
ionized  uranium. 

An  accurate  theoretical  model  of  the  uranium  atom  would  represent 
a  monumental  step  in  understanding  complex  atomic  systems.  The  mathema- 
tical description  of  many-electron  atoms  has  been  attempted  using  several 
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methods  imbedded  with  various  approximations.      In  many  situations, 

experimental  results  are  available  to  substantiate  theory  or  at  least  to 
raise  questions  about  the  validity  of  certain  aspects  of  theoretical 
treatment.  Currently,  there  is  little  conclusive  theoretical-experimental 
validation  of  the  atomic  properties  of  uranium  and  uranium  plasmas.  How- 
ever, there  is  a  significant  effort  under  way  by  Steinhaus  et  al.   ' 

to  establish  experimentally  the  energy  levels  of  neutral  and  once- 

20 
ionized  (UI  and  DII)  uranium.  This  effort  evolved  from  work  by  Schuunnans 

and  Kiess  et  al."'  in  1946.  With  subsequent  improvement  in  opti cs , spectro- 
scopic techniques,  and  atomic  structure  calculations,  confidence  is  increas- 
ing in  energy-level  definition  but  progress  is  slow. 

The  specific  emission  coefficient  e;>^(T,P)  characterizes  the  light 
emission  from  a  volume  element  of  plasma.  It  is  defined  as  the  amount 
of  energy  emitted  per  unit  time,  from  a  given  volume,  into  a  specified 
solid  angle  and  wavelength  interval.  It  can  be  composed  of  continuous 
and/or  discrete  components,  both  of  which  are  strong  functions  of  plasma 
temperature  and  density  (pressure)  of  the  radiating  species.  The  form- 
ulation of  an  adequate  theoretical  model  of  the  uranium  atom  can  be 
strongly  assisted  (perhaps  out  of  necessity)  by  obtaining  from  experi- 
ment detailed  knowledge  of  the  uranium  emission  coefficient.  Because 
calculation  of  such  a  property  (for  uranium)  is  impossible  without  a 
model  that  suffers  from  several  approximations,  experimental  verfication 
through  emission  coefficient  measurement  is  needed.  This  thesis  reports 
on  an  experiment  designed  to  measure  the  specific  emission  coefficient 
of  a  uranium  plasma  and  to  relate  these  data  to  state-of-the-art 
theoretical  predictions. 
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Miller   and  Marteney  et  al .   have  acquired  emission  coefficient 
data  from  UF^  shock-tube  and  radio-frequency  induction-heated  Ar-UF, 
plasmas,  respectively.  Uranium  plasma  emission  coefficients  obtained 
from  UF^  discharges  present  distinct  impurity  problems,  potentially 
resulting  in  distorted  emission  coefficient  wavelength  dependence. 
A  dc  uranium  arc  was  chosen  as  the  light  source  for  the  present  ex- 
periment to  r-educe  plasma  impurities  and  provide  a  steady-state,  less 
contaminated  plasma  to  determine  its  emission  properties.  The  wave- 

o  o 

length  bandwidth  considered  is  from  1050  A  to  an  upper  limit  of  6000  A. 
The  emission  data  are  then  compared  to  similar  data  generated  by  other 
research  groups.  Also  included  is  a  comparison  to  theoretical  emission 
coefficient  predictions.  In  summary,  this  effort  was  conceived  to  estab- 
lish a  unified  picture  of  the  progress  made  over  the  past  few  years 
concerning  the  experimental  and  theoretical  investigations  characterizing 
uranium  plasma  emission. 

Chapter  II  describes  the  plasma  diagnostics  necessary  for  the 
determination  of  basic  plasma  properties  such  as  temperature,  particle 
densities,  and  radiation.  Particular  attention  is  given  to  the  applica- 
tion of  such  diagnostics  to  uranium  plasmas.  Uranium  arc  plasma  genera- 
tion and  arc  stability  are  discussed  in  Chapter  III.  Also  examined  are 
the  various  methods  of  data  acquisition  and  intensity  calibration  applied 
in  the  course  of  this  investigation.  Chapters  IV  and  V  indicate  the 
emission  measurements  for  uranium  plasmas  at  two  distinct  arc  conditions. 
Brief  descriptions  and  comparisons  of  other  similar  experiments  are 
given  amoung  these,  the  present  effort,  and  theoretical  predictions. 
In  Chapter  VI  theoretical  models  of  the  uranium  atom  are  critically 
considered.  A  comparison  is  made  between  emission  peak  locations  of  the 


present  results  to  theoretically  predicted  locations  for  the  higher 
temperature  arc  plasma.  A  discussion  of  the  major  points  of  this  over- 
all effort  then  concludes  this  study. 


IT.      'PT.ASm   TURORY  AND  DIAGNOSTICS 

Plasma  temperature  and  density  diagnostics  are  based  on  the 
assumption  that  relevant  information  can  be  extracted  from  theoretical 
descriptions  of  microscopic  plasma  processes.  Generally,  plasma 
constituents  exhibit  balance  between  the  population  and  depopulation 
of  neutral  and/or  ionic  energy  levels.  This  implies  that  for  every 
excitation  event  (to  use  an  example)  there  is  a  corresponding  de- 
excitation,  not  necessarily  brought  about  by  the  inverse  excitation 
process.  Plasma  particles  that  are  excited  and  de-excited  by  the 
same  mechanism  demonstrate  detailed  balancing.  Arc  plasmas  are  usually 
considered  to  be  collision-dominated  in  the  sense  that  e  -atom  colli- 
sions cause  most  of  the  excitations  and  de-excitations.  However, 
radiative  de-excitation  can  also  be  important  as  indicated  by  strong 
photon  output.  Therefore,  arc  plasmas  are  rarely  characterized  by  de- 
tailed balancing,  and  various  approximations  such  as  local  thermo- 

24 
dynamic  equilibrium  (LTE)  and  partical  LTE   must  be  used. 

II- 1 .      Equi  I i b Y'i wn_  Considerations 

Griem^^  has  developed  criteria  that  indicate  which  equilibrium 

state  applies  to  a  given  plasma.  Most  of  these  tests  set  limits  on 

the  dominance  (or  lack  thereof)  of  collisional  rates  over  radiative 

rates.  Griem  developed  the  following  criterion  for  complete  LTE  in 

a  hydrogen- like  system: 


e   -    -   »  .Tr-|  1  .2 


Z^H  /  \  ^  ^H 


_3 
Where    n  =  electron  density  (cm  ), 

Z  =  net  plasma  charge  (Z=l  for  neutral;  Z=2 
for  singly  ionized) , 
/^[-     =  energy  (eV)  of  1st  excited  state  with  respect 
to  the  ground  state, 
kT    =  thermal  energy  of  plasma  electrons  (eV),  and 
E,,   =  ionization  potential  of  hydrogen  (eV). 
This  criterion  estimates  the  lower  limit  of  electron  density  re- 
quired to  maintain  a  hydrogenic  system  in  complete  LTE.  A  hydrogenic 
system  is  usually  characterized  by  a  large  energy  gap  between  the 
ground  state  and  first  excited  level,  and  successively  smaller  energy 
separations  between  levels  as  continuum  is  approached.  Hydrogen  (and 
helium)  is  also  one  of  the  more  difficult  elements  to  bring  into  LTE 
because  of  its  relatively  high-lying  first  excitation  level.  Thus, 
hydrogen  can  usually  serve  as  an  upper-1 imit  estimation  for  the  validity 
of  LTE  for  more  complex  systems.  However,  the  electron  density 
criterion  offered  by  Eq.  (II-l)  is  not  directly  applicable  to  the 
uranium  system  because  the  known  energy-level  description  of  the 
uranium  atom  (or  ion)  simply  does  not  fit  the  hydrogenic  picture.  In 
fact,  the  first  excitation  levels  of  UI  and  UII  are  ^  0.077  eV  (UI) 
and  0.03  eV  (UII),  which,  for  singly  ionized  uranium  at  8500  K,  implies 

n  -  3  9  X  10^  cnf"^.  This  is  well  below  any  reasonable  values  for  n 
e  ^ 

basically  because  of  the  strong  dependence  on  AEq^  which  has  a  very 
small  value  for  neutral  and  once-ionized  uranium.  Applied  to  hydrogen 

at  8500  K,  n  ^  8.8  x  10^^  cm"^  this  limiting  value  for  n^  is  probably 
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too  restrictive  for  the  uranium  system.  Griem,   McWhirter,   and 


Wilson^'^  have  developed  a  test  which  depends  on  the  energy  level 
structure  of  a  specific  atomic  system,  thus  reducing  the  necessity 
of  assuming  hydrogen-like  characteristics.  This  criterion  (hereafter 
called  the  ladder  criterion)  is  more  applicable  to  complex  systems. 
The  ladder  criterion  presumes  that  the  most  difficult  level  to 
populate  will  form  the  largest  energy  gap  ascending  the  energy  level 
diagram.  Existence  of  LTE  in  an  optically  thin  collision-dominated 
plasma  is  tested  by 

n^>  C  T^^  (E^-  E.)3  ,  (II-2) 

where  C  =  1.6  x  10^"^  (McWhirter^'') . 
The  electron  temperature  T^  has  units  of  K;  the  largest  energy  gap  in 
eV  for  the  atom  under  study  is  (E^^  -  E^. ).  Equation  (II-2)  is  not  suf- 
ficient for  LTE  because  e"  relaxation  time,  e"-atom  relaxation  time, 
and  atom  relaxation  time  must  dominate  whatever  other  transient  phe- 
nomena occur  in  the  plasma.  Most  stable  steady-state  arcs  satisfy 
these  relaxation  requirements. 

To  apply  Eq.  (II-2)  to  a  uranium  plasma  it  is  necessary  to  define 
the  largest  energy-level  gap  found  in  its  level  diagram.  Complete 
tabulation  of  this  information  is  unavailable,  but  a  representative 
gap  width  may  be  defined  by  examination  of  the  work  by  Steinhaus 
et  al.^*^'^^  and  Blaise  and  Radziemski.    Considering  neutral  uranium 
(UI),  the  620.323  cm"^  ->  3800.829  cm"^  levels  provide  the  largest  gap 
('■^  3180  cm"^).  This  implies  an  electron  density  of  7.2  x  10   cm 
at  5500  K.  Similarly  for  singly  ionized  uranium  (UII),  the  largest 
gap  is  about  2126  aif^  (2294.70  ->   4420.87).^°  Assuming  an  electron 
temperature  of  10  000  K  (UII),  the  corresponding  electron  density 


limit  is  then  2.93  x  10   cm"  .   If  partial  LTE  is  assumed,  the  larg- 
est energy  gap  need  only  be  defined  above  the  thermal  energy  limit. 

In  conclusion,  although  there  is  no  feasible  way  to  determine 
accurately  the  LTE  limits  on  electron  density  for  atomic  and  once- 
ionized  uranium,  reasonable  estimates  can  be  obtained  using  Eqs.  (II--1)- 
and  (II-2). 
17-  2 ._ Plas ma  Radiation 

Many  factors  affect  the  net  photon  output  from  a  plasma.  Optical 
depth  and  absorption  are  the  most  important.   In  many  laboratory 
plasmas,  photons  emitted  from  the  plasma  central  region  may  be  absorbed 
en  route  to  the  outer  boundary.  The  degree  of  trapping  implies  an 
optical  depth,  t,^(x),  which  is  simply  the  effective  absorption  coeffi- 
cient  integral  k'(x)  over  a  given  homogeneous  plasma  depth  X.   The 
relation  between  the  intensity  I^(x)  and  energy  emitted  as  a  function 
of  the  line-of-sight  depth  into  the  plasma  is  governed  by  Eq.  (II-3), 
the  equation  of  radiative  transfer. 

dl  (x)         z   (x) 

-X--v=  I  (X)  ---^-rvT  •  (^^-3) 

dr  (x)    \)  K  (xj 

The  specific  emission  coefficient,  e,^(x),  (energy/volume-wavelength- 
time-solid  angle)  and  c  (x)/k:^(x)  is  the  source  function  of  radiation 
(energy/area-wavelength-time-solid  angle).  The  source  function  for  an 
LTE  plasma  is  equivalent  to  the  Planck  function. 


In  this  sense  X  connotes  a  homogeneous  plasma  depth,  and  the 
subscript  V)  implies  that  there  is  a  spectral  dependence  to 

T  and  ^:. 


The  general  solution  to  the  equation  of  radiative  transfer  is 

X 

I  (x)  =  I  (o)  e-"v^^^  +  f     B  (T)  k;(x')  e'^^^^^  dx^  ,  (II-4) 

Q 

where  I  (x)  =  observed  intensity  at  X,  v  (energy/area-time-solid 
angle-wavelength, 
I  (o)  =  intensity  of  background  source  (if  there  is  one)  at  v, 
B  (T)  =  Planck  function  at  T,  v  (energy/area-time-solid 
angle-wavelength) , 

K-'U')  =   effective  absorption  coefficient  (length"  )  at  x^,  v, 
and 

e   (x')  =  specific  emission  coefficient  at  x^  v  (energy/volume- 
tinie-solid  angle-wavelength). 

In  Eq.  (II-4),  a  situation  is  possible  where  aJl  emitted  photons 
transport  beyond  the  outer  plasma  boundary,  that  is,  if  t^(x)  ^  o. 
Then  the  observed  intensity  is  simply  the  emission  coefficient  integral 
over  the  1 ine-of-sight  depth,  provided  the  plasma  is  homogeneous. 
This  plasma  is  designated  optically  thin.  The  other  extreme  would  be 
T  (x)  -y  «s  here  all  emitted  photons  are  trapped  within  the  outer  plasma 
boundary,  and  then  the  plasma  radiates  as  a  surface  where  I^  ->  B,^. 
This  plasma  is  designated  optically  thick  and  exhibits  a  blackbody 
spectral  distribution  of  radiation. 

Generally,  in  relating  the  concept  of  optical  depth  to  arc  plas- 
mas, the  line-of-sight  depth  and  the  effective  absorption  coefficient 
must  be  determined.  Usually  the  line-of-sight  depth  will  be  quite 
small  (for  arcs),  making  the  effective  absorption  coefficient  the 
important  quantity.  Unless  the  plasma  pressure  is  50  atmospheres  or 
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greater  (dense  arc  plasmas),  the  only  significant  absorption  will  be 
located  (on  the  wavelength  scale)  at  some  of  the  spectral  line  cen- 
ters, usually  formed  in  the  lower  energy  levels.  In  many  cases  for 
low-  and  intermediate-density  plasmas,  the  effective  absorption 
coefficient  will  be  small  for  most  wavelengths. 

Line,  recombination,  and  continuum  radiation  are  three  basic 
radiation  types  occurring  in  an  arc  plasma.  Line  radiation  is  usually 
associated  with  relatively  low  temperatures,  i.e. ,4000  K-15  000  K; 
recombination  and  continuum  radiation  can  be  substantial  when  the 

characteristic  plasma  temperature  is  >  15  000  K.  Detailed  treatments 

31 
of  recombination  and  continuum  radiation  are  given  by  Cooper. 

Spontaneous  emission  and  stimulated  emission  result  in  line  radi- 
ation. Stimulated  emission  (often  thought  of  as  negative  absorption) 
is  a  difficult  item  to  isolate  and  therefore,  is  usually  defined  as 
an  effective  absorption  coefficient  for  any  plasma  as  shown  by: 

K'(x)  =  k:   (x)  +  K^   (x)  (II-5a) 


V 


(x)  =  K^     (x)  -  Kg^  (x)    , 


where  k"  =  total  (measured)  effective  absorption  coefficient  at 

-1 

x ,  V  cm   , 

K,"  (x)  =  total  effective  line  absorption  coefficient  at  x,  v, 

V 

K,  (x)  =  line  absorption  coefficient  at  x,  v, 

V 

K  (x)  =  total  continuum  absorption  coefficient  at  x,  v,  and 
c 

V 
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K   (x)  =  negative  absorption  coefficient  caused  by  stimulated 


emission  at  x,  v. 
Thus,  Kc-r  (x)  represents  an  effective  decrease  in  k'(x).  If  an  energy 

V 

absorption  transition  by  a  bound  electron  of  an  atom  is  considered, 
where  u  =•>  upper  level  and  £  =>  lower  level,  the  description  of  line 
absorption  at  a  given  frequency  is  given  by 


hv 


K^  (x) 

V 


l-^u 


"n  B„     *  (v) 


(II-6) 


where 


Vn       =  the  frequency  of  the  transition   (s     ), 

_3 
n     =  population  density  of  the  lower  state   (cm     ): 

B,       =  Einstein's  probability   for  absorption   (cm/g) 

(|)Jv)   =  absorption  line-shape  function  with 

d)    (v)   dv  =   1  ,   and 

K 

me 


i 


X  =  line-of-sight  distance  into  plasma. 
Similarly,  the  line  emission  coefficient  for  a  homogeneous  plasma 
at  frequency  v  is  described  by 

hv, 


e^  (x) 


u-^y 


4'^ 


u   u    ^£^  ' 


II-7) 


where  £.,  (x)  is  the  specific  line  emission  coefficient  at  v  and  plasma 

V 

depth  x.  Actually,  a  continuum  emission  coefficient  should  also  be 
considered  for  completeness,  but  this  is  not  included  here  because  of 
its  small  (or  inextractable)  contribution.  The  population  density 
(cm  )  of  the  spontaneous  upper  level  transition  and  Einstein  probabil- 
ity for  such  a  transition  are  given  by  n^  and  A^_^^(s  ),  respectively. 
The  normalized  line-shape  function  is  indicated  by  cf>  (v)  ,  such  that 
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A.  (b   (v)dv  =  1.  The  line  intensity  I  (x)  from  such  a  transition 
line  "^e.  '  V 

v;ould  simply  be  determined  by  multiplying  the  line  emission  coefficient 
by  the  appropriate  LOS  plasma  depth  through  a  constant  emission  zone. 

Arc  plasmas  are  inherently  inhomogeneous  to  varying  degrees.  Free- 
burning  arcs  usually  exist  with  significant  temperature  gradients  along 
the  major  part  of  an  approximately  cylindrical  arc  column  radius,  which 
results  in  radially  varying  intensities.  Therefore,  the  emission  coef- 
ficient will  also  have  a  radial  dependence  that  must  be  extracted  by 
unfolding  methods  applied  to  observed  intensities.  Many  wall-stabilized 
arcs  have  less  severe  temperature  gradients,  with  the  exception  of  that 
region  approaching  the  wall."^^"'^^  This  implies  the  possibility  of  an 
approximately  homogeneous  nature  in  temperature  and  density  in  the  major 
part  of  the  cylindrical  arc  plasma  in  the  radial  direction.  Although 
a  relatively  constant  temperature  profile  does  not  assure  homogeneous 
density  and  emission  profiles,  it  can  be  a  strong  indicator. 
jT-g.  mission  Coefficient  Detei'mination 

The  geometry  of  inhomogeneous  optically  thin  arc  plasmas  is  approxi- 
mated by  a  number  of  concentric  zones  about  the  vertical  axis,  as  shown 
in  Fig.  IT-l.  Each  zone  is  assumed  to  display  constant  emission  which 
is  a  function  of  a  single  temperature  and  density  for  that  zone.  The 
intensity  at  a  given  location  along  the  arc  chord  is  the  sum  of  emission 
contributions  for  each  zone  intercepted  by  a  LOS  ray  passing  through  the 
geometry.  The  emission  coefficient  must  then  be  unfolded  from  the 

measured  integrated  intensities  by  the  familiar  Abel  transform.  The 

38  . 
analytical  form  of  the  transform  equations   is 
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i(y)  =  2 


X 

f    e(r)dx  =  2  / 


r 

^   e(r)rdr  ,  and         (II-8a) 


V' 


.,  2  2 
y  "Vr  -y 


e(r)  =  -i  rriy)^ 

'""  7   J  2  2 
r  ^y  -r 


;il-8b) 


The  geometrical  relations  are  developed  using  Fig.  II-l. 

There  are  two  common  approaches  used  in  solving  Eqs.  (II-8a)  and 
(II-8b)  for  e(r).  The  first  is  by  fitting  an  appropriate  polynomial 
to  the  experimental  intensities,  differentiating,  and  using  Eq.  (II-8b) 
to  obtain  e(r).  The  details  of  such  a  method  are  given  in  Ref.  39. 
The  second  method  approximates  the  integrals  with  sums  and  extracts  the 
desired  information  numerically.  The  numerical  form  of  Eq.  (II-8a)  is 
given  by 

I^(y)  =  2  y  f.  .(r)  £.  .   .  (11-9) 

J 

In  Fig.  II-l,  the  i^^  LOS  in  the  j^*^  concentric  ring  is  defined  by  the 
length  coefficients,  I..,   and  e-(r)  represents  the  average  emission 
coefficient  for  the  j^^  ring.  Equation  (II-9)  results  in  a  system  of 
simultaneous  linear  equations  for  each  c.(r),  I-(y),  which  can  be 
solved  by  matrix  inversion  to  obtain  the  e.(r)  vector. 

When  the  plasma  under  investigation  has  significant  intensity  gra- 
dients, a  number  of  rings  (as  many  as  40)  may  be  necessary  to  approximate 
adequately  the  intensity  profile  and  allow  computation  of  emission 
coefficients  within  acceptable  error.  Each  ring  requires  a  correspond- 
ing intensity  determination,  which  may  be  too  difficult  to  obtain, 
depending  on  plasma  stability.  However,  too  many  subdivisions  may  cause 
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oscillations  in  the  unfolding  process  if  an  exceptionally  smooth  in- 
tensity profile  is  not  available.  Substantial  error  (20  -->  30%)  on  the 
intensity  profile  measurement  (especially  for  the  outermost  ring)  v\;ill 
be  propagated  to  the  central  zone  and  will  yield  a  poor  estimate  of  its 
emission  value.  For  this  situation,  too  many  rings  will  introduce  large 
oscillations  and  an  unacceptable  error  propagation.  Obviously,  there  is 
an  optimum  number  of  rings  depending  upon  the  overall  intensity  profile 
shape  and  the  experimental  error  value,  particularly  at  the  outermost 
ring. 

Many  studies  have  dealt  with  various  techniques  to  implement  the 

inversion  "   but  these  primarily  address  algorithmic  problems  (inter- 

44       42  43 

polation,  smoothings,  etc).  Kock,   Nestor,   and  Bockasten 

address  error  propagation  by  using  numerical  methods  for  inversion. 

These  analyses  imply  that  the  correct  number  of  subdivisions  is  a  direct 

function  of  the  experimental  error  in  the  outermost  zone  as  well  as 

44 
the  shape  of  the  intensity  profile.    When  error  in  the  outer  zones 

approaches  20%  or  more,  four  or  five  zones  may  define  the  limiting 
accuracy. 

To  relate  these  considerations  to  the  inversion  of  the  wall- 
stabilized  uranium  arc  intensity  data,  a  four-ring  numerical  unfolding 
scheme  was  chosen  because  of  time  limitations  involved  in  acquiring 
zonal  intensity  data  for  more  than  four  zones,  and  the  experimental  er- 
ror associated  with  the  outermost  zonal  intensity  was  -^  +25%.  When  this 
error  is  considered  with  a  10%  calibration  error  and  an  estimate  of 
error  for  using  only  four  rings,  the  total  error  estimate  of  central 
zone  emission  coefficients  ranges  from  28  -^  36%. 


II- 4.      Plasma  Temperature 
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Many  techniques  used  for  temperature  diagnostics  were  developed  by 
Griem,^"^  Hefferlin,^^  and  Lochte-Hol tgreven.    The  following  discussion 
considers  only  those  temperature  diagnostics  applied  to  a  uranium  arc 
plasma,  including  the  Boltzmann  plot,  norm  temperature,  relative  norm 
temperature,  and  the  modified  brightness-emissivity  methods. 

A. Boltsmann  Plot  Method 

•k 

The  integrated  line  intensity  for  a  homogeneous,  optically  thin 
plasma  into  a  depth  x  is 

hv 


I 


U'^l 


A  „  n  X 


\)    4tt     u->£.  u 


:ii-io) 


If  the  plasma  is  in  complete  LTE  and  inhomogeneous,  n^(T)  is  determined 
by  the  Boltzmann  factor  yielding 


I     = 


hv     „     g 

u->s;.       u 


47T 


zm    V^%(T)   e       u.       £.. 


11-11) 


where       T 
A 


excitation  temperature   (K)    , 


•1 


V 


u->-ii 


n,(T) 


Z(T) 

E 
u 

I. . 


transition  propability  (seconds  )  , 

transition  frequency  , 

statistical  weight  of  the  upper  level  , 

ground  state  population  density  of  a  particular  ionization 

state  (#/cm^-)  , 

partition  function  , 

energy  of  upper  level  (eV),  and 

ring  depth  of  zone  of  assumed  constant  emission  (Fig.  II-l), 


'integrated  in  this  sense  means  integration  over  the  line-shape  function 


^ 


(|)  (v)dv  has  already  been  performed, 
line  t: 
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Taking  the  natural  log  of  both  sides  of  Eq.  (11-11)  and  rearrang- 
ing terms  yield  an  equation  analogous  to  that  of  a  straight  line: 


'  „  I  )      (h  n  (T)ii.  .  i    (  T 


Hrei^N--^Ty-r  i-r--     '"-"= 


If  for  tv/o  or  more  spectral  lines,  In  (constant  *  I  )  is  plotted    . 

against  E   ("m"  refers  to  a  particular  transition),  the  resultant 
^      urn  ^  ^ 

slope  of  the  curve  will  approximate  -1/kT.   If  I  is  measured  in  abso- 
lute units,  the  ordinate  intercept  also  determines  the  LIE  ground  state 
number  density  for  that  species. 

Equation  (11-12)  was  developed  considering  complete  LTE  and  can  be 
modified  for  use  with  less  restrictive  equilibrium  concepts.  When 
partial  LTE  is  assumed,  the  Boltzmann  plot  method  for  complete  LTE  is 
modified  to  yield 

where:  E  =  energy  level  of  lower  energy,  and 

g  =  lower  level  statistical  weight. 
The  essential  difference  between  Eq.  (11-12)  and  Eq.  (11-13)  is  that 
the  intercept  term  no  longer  defines  the  true  ground-state  density  in 
partial  LTE. 

Application  of  the  Boltzmann  plot  technique  for  relative  or  abso- 
lute temperature  determination  requires  the  following  conditions: 

(1)  at  least  partial  LTE  must  exist  in  an  optically  thin  plasma, 

(2)  many  transitions  should  be  used. 
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(3)  these  spectral  lines  should  have  a  correspondingly  significant 
energy-level  separation,  and 

(4)  this  technique  can  be  used  for  inhoinogeneous  plasmas,  but  the 
Intensity  data  ideally  should  first  be  spatially  resolved  to 
express  the  emission  coefficient  as  a  function  of  plasma  radius 
(in  cylindrical  geometry).  If  there  is  no  spatial  resolution 
the  computed  temperature  will  be  an  approximate  value,  but  in 
many  instances  it  will  yield  a  good  estimate  of  an  average  arc 
temperature.  The  accuracy  of  these  temperatures  is  governed  by 
intensity  measurement  and  uncertainty  in  atomic  constants.  The 
largest  source  of  error  for  uranium  plasmas  is  usually  caused 
by  the  uncertainty  of  gA  values  which  can  be  as  much  as  50%. 
Only  in  extreme  cases  will  the  average  centerline  temperature 
deviate  by  more  than  20%  from  the  true  centerline  temperature. 

To  apply  this  method  to  uranium  plasmas  (or  any  method  dependent  on 
gA  values),  availability  of  well-defined  spectral  lines  is  severely 
limited.  These  lines  must  meet  the  previously  mentioned  requirements 
and  be  locatable  with  the  instrumentation  to  be  used.  Application  of 
the  Boltzmann  plot  method  to  uranium  plasmas  is  straightforward  in 
principle;  but  practically  speaking,  a  very   tedious  task  strongly  de- 
pendent on  spectrograph  resolution,  line  identification,  and  availa- 
bility of  relevant  constants. 

B.        Nor'm-Temperatui'e  Method 

If  the  equilibrium  emission  e,  (T)  of  a  transition  in  a  specified 
ionic  state  is  plotted  against  temperature,  and  if  the  temperature  is 
high  enough,  the  resultant  curve  will  be  peaked  at  the  norm-temperature 
T  ,  defined  on  an  eniission-vs-temperature  plot  as  the  temperature 
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at  which  d£  (T)/dT  =  £"(T)  =  0.  The  norm  temperature  can  then  be  used 
to  estimate  the  characteristic  plasma  temperature.   First,  it  is  neces- 
sary to  determine  which  ionization  state  corresponds  to  identified 
emission  'lines  generated  by  the  plasma.  Then  e  (T)  is  computed  for  a 
spectral  line  known  to  be  from  the  dominant  ionic  stage  by  using  appro- 
priate equilibrium  relations.  Then  e'(T)  =  0  is  determined,  thus  de- 
fining T  .  The  maximum  temperature  of  an  equilibrium  plasma  emitting 
a  line  in  a  particular  dominant  ionic  state  can  be  on  the  order  of  T  , 
which  is  indicative  of  the  energy  necessary  to  ionize  the  plasma  to  that 
state. 

Obvious  disadvantages  of  such  a  method  are  that  T  is  only  an 
estimate  of  the  characteristic  plasma  temperature,  and  while  transition 
probabilities  need  not  be  known,  the  Saha  equation  must  be  solved  to 
obtain  the  neutral  and  ionic  number  densities,  thereby  requiring  LTE. 
The  priinar^y  advantages  are  that  only  spectral  line  is  necessary  and  that 
this  method  can  easily  be  adapted  to  the  relative  norm-temperature  method. 

^1  Ralative  Temperature  Method 

The  relative  temperature  method  is  an  extension  (in  many  cases)  of 
the  norm- temperature  method  in  that  the  norm-temperature  can  be  used  as 
a  reference  plasma  temperature  and  temperatures  at  other  spatial  loca- 
tions related  to  it.   It  can  be  used  for  plasmas  which  have  cylindrical 
geometry  with  radial  temperature  profiles  T(r)  and  at  least  one  defined 
temperature  such  as  T  .   If  this  is  the  case,  the  resultant  emission 
profile  c  (r)  pertaining  to  a  specific  transition  can  be  similar  to  that 
shown  in  Fig.  11-2;  one  point  on  the  desired  temperature  profile  T(r) 

is  defined  by  c  (r  )  and  T  .   From  equilibrium  relations,  ratios  of 
■^  V  0^     n        ^ 

e,(r  )/c  (r.)  may  be  computed  and  the  corresponding  T(r.)  determined. 
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e  (r.)   Z[T(r  )]  n^(r.) 


S^'o^   Z[T(rJ]  nT(r  ) 


exp 


U 

k 

L      J 

J(r.) 


1 

T(r 


0' 


(11-14) 


Any  temperature  on  T(r),  determined  by  other  independent  methods, 

will  suffice  for  T(r  ). 

^  o' 

The  relative  temperature  method  is  particularly  straightforward 
wlien  dealing  with  the  norm-temperature  situation.  Transition  probabili- 
ties are  not  needed,  but  the  temperature  dependence  of  partition  function 
should  be  considered.  Only  one  transition  need  be  used,  thus  eliminating 
any  calibration  procedures.  This  method  is  very  attractive  for  use  with 
spectroscopical ly  complex  plasmas  because  it  is  independent  of  gA  values. 

!!li^—^H2ttklji: et^^  Bi'-ightness-Emissivitii  Method 

Temperature  measurement  methods  previously  discussed  require  the 
assumption  of  negligible  absorption  or  an  optically  thin  plasma.   If 
certain  spectral  lines  emitted  by  a  thin  plasma  have  measurable  absorp- 
tion (say  at  the  spectral -line  centers),  temperature  may  be  determined 
using  spectral-line  absorption  as  a  basis. 

The  brightness-emissi vi ty  method  (BEM)  adapts  the  radiative- 
transfer  equation  (Eq.  (II-3)]  solution  to  a  homogeneous  LTE  absorbing 
plasma.'  '    The  following  equations  form  the  basis  for  this  method. 


I'(x)  =  I"(x)  e"''^^^''^  +  I^(x) 


V 


I^(x) 


i^(x: 

V 


1^(0)  e-^v(x)  ,  bP(,; 


-T,,(X) 


11-15' 


1^,(0)  e  ^^ 


r.,(x)    p 

+  B  T  e 


Emissivity,  as  considered  in  this  instance,  is  the  ratio  of  the  plasma 
intensity  to  the  Planck  intensity  at  a  specified  temperature  and  wave- 
length. 
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where     'u  =  emissivity, 

X  =  LOS  plasma  depth, 

T 

I    (x)   =  total   observed  intensity  at   (v,   x), 

\) 

I^{x)  =  total  external  source  intensity  without  absorption 

V 

at  (v,  x), 
1^(0)  =  total  external  source  intensity  without  absorption 

at  V),  x  =  0  (the  plasma  boundary),  and 
B^(T)  =  Planck  function  indicative  of  the  brightness  temperature 
of  the  plasma. 
The  known  quantities  I."^(x),  Iy(x),  'r^(x),  and  e^  allow  determina- 
tion of  the  plasma  brightness  temperature  using  the  Planck  function 
for  the  plasma. 

The  primary  disadvantage  of  the  BEM  is  lack  of  applicability  to 
inhomogeneous  plasmas  such  as  those  formed  with  arcs.  Usher  and 
Campbell "^^  have  adapted  the  BEM  to  the  homogeneous  case.  The  modified 
BEM  differs  from  the  BEM  in  that  an  unfolding  scheme  resolves  the  absorp- 
tion coefficient  and  temperature  profiles  spatially. 

The  absorption  coefficient  of  a  spectral  line  is  determined  by  a 
constant  intensity  background  source  and  measurement  of  the  wavelength 
attenuation  as  the  line  passes  through  the  plasma.  If  the  inhomogeneous 
plasma  is  composed  of  homogeneous  rings  (similar  to  Fig.  II-l),  an  un- 
folding technique'^^  can  be  used  to  calculate  the  average  line-center 
absorption  coefficient,  k1{X),   for  each  ring  from 
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T       P 
Voltage  signals  V^,  V',  and  V.  are  from  a  photoniul tiplier  tube  for  the 

i^'^  LOS  position.  The  voltages  indicate  background  source  intensity, 
total  attenuated  intensity,  and  plasma  intensity,  in  that  order.  The 
£..  in  Eq.  (11-16)  represents  the  length  segments  in  the  j   ring  along 
the  i^^^  LOS  position  as  shown  in  Fig.  II-l.   Figure  II-3  shows  a  typical 
oscillogram  of  the  photoniul tipl ier  output  related  to  Eq.  (11-16).  The 
background  source  of  known  characteristics  (Xenon  flashtube)  is  flashed 
on  the  line-center  of  interest  as  shown  on  the  upper  trace.  The  lower 
trace  shows  the  flashtube  signal  spread  in  time  to  facilitate  intensity 
voltage  measurement.  The  oscillograms  are  recorded  at  different  chordal 
positions  of  the  plasma  and  Eq.  (11-16)  is  used  to  determine  the  line- 
center  absorption  profile  across  the  radial  dimension  of  an  assumed 
cylindrically  symmetrical  arc. 

The  temperature  is  determined  by  an  extension  of  the  brightness- 
emissivity  method  to  the  inhomogeneous  case.  The  technique  uses  measured 
voltages  and  calculated  line-absorption  values.  The  temperature  method 
requires  that  the  background  source  temperature  (or  equivalently,  the 
intensity)  be  known.  The  average  temperature  in  the  j   ring  of  the 

plasma  is  determined  from 
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The  Planck  functions  B<.  and  B-  represent  intensities  of  the  background 
source  and  the  j*"^^  plasma  ring  at  wavelength  X,   and  Tj^  is  the  brightness 
temperature  of  the  background  source.  The  B.'s.are  calculated  using  the 

J 


me 


asured  voltages  and  the  computed  plasma  absorption  profile.  The 


complicated  expression  for  B-  is  found  in  Ref.  47. 
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There  is  very  little  available  information  about  atomic  properties 
of  uranium.  This  information  would  be  useful  when  performing  plasma 
temperature  diagnostics  on  spectroscopically  complex  elements;  however, 
the  modified  BEM  is  particularly  effective  for  use  with  uranium  plasmas 
because  it  is  independent  of  these  constants.  Because  this  technique 
deals  adequately  with  plasma  inhomogeneities,  it  is  suitable  for  appli- 
cation to  arc  plasmas.  Implementation  of  the  method  is  cumbersome 
because  it  means  a  long-duration  ('x.  minutes)  steady-state  plasma,  a 
well-defined  standard  background  source  and  its  associated  circuitry  to 
provide  a  rather  elaborate  data  acquisition  sequence. 
I I- 5.      Dens i ty  Msasuvements 

To  define  plasma  pressure,  total  or  partial,  it  is  necessary  to 
determine  ground-state  population  densities  for  all  ionization  stages 
in  the  plasma.  The  equation  of  state  (in  conjunction  with  Dalton's  law 
of  partial  pressure)  for  each  plasma  component  is  then  used  to  define 
a  total  pressure.  A  substantial  effort  is  involved  in  establishing 
population  densities.  The  determinations  must  rely  on  an  accurate  atomic 
description  of  the  plasma  constituents  and  precise  intensity  calibration. 
Griem^^  states  that  even  in  an  optimum  situation,  it  is  often  impossible 
to  reduce  density  error  below  30%.  Therefore,  most  plasma  density 
determinations  are  order-of-magnitude  estimates. 

Two  density  diagnostics,  which  were  applied  to  the  uranium  arc 
plasma,  are  described  in  the  following  sections:  absolute  line  inten- 
sity and  pressure-temperature  correlation  (PTC). 

/L Absolute  Line  Intensity  Method 

The  integrated  line  intensity  for  a  transition  from  level  u  to 
level  £  is  given  by  Eq.  (11-10).  Equation  (11-10)  applies  to  a  line 
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emitted  tYom  a  honiegeneous  optically  thin  LTE  plasma.  Modification  to 
an  absorbing  plasma  requires  compensation  (build-up)  along  the  entire 
line-shape  profile.  Obtaining  absorption  build-up  factors  at  many 
points  along  the  line  profile  may  not  be  possible  because  of  line-wing 
overlap,  especially  in  plasmas  displaying  complex  spectra  where 
isolated  lines  may  not  exist.  Fortunately,  absorption  in  many  arc 
plasmas  tends  to  be  concenty^ated  at  spectral -1  ine  centers  (on  the  wave- 
length scale).  Therefore,  an  absorption  coefficient  determined  at  the 
spectral-line  center  approximates,  to  a  usually  acceptable  degree,  the 
maximum  absorption  coefficient  taken  over  the  entire  line  profile.  This 
remains  a  good  approximation  as  long  as  the  spectral  line  exhibits  a 
sharp  profile.  Using  the  above  as  a  basis,  absorption  build-up  may  be 
incorporated  into  the  line-intensity  equation  by  the  following: 

hv  ,.       G         -E  /kT^ 
l'     =  -^      -^-h         n  e  '   ^  I,.   Y  ,       (11-18) 
^total    ^'"     Z   u    0         13 

where  y  represents  a  dimensionless  absorption  build-up  factor. 

The  variables  in  Eq.  (11-18)  are  the  intensity,  the  excitation 
temperature,  and  the  ground-state  number  density.  Clearly,  to  find  n^ 
for  a  homogeneous  plasma,  absolute  intensity  units  must  be  known  and 
excitation  temperature  must  be  determined  by  an  independent  method.  For 
an  inhomogeneous  plasma,  unfolding  should  be  performed,  or  the  calculated 
number  density  would  indicate  an  approximate  value.  The  absorption 
build-up  is  determined  experimentally,  but  this  factor  will  be  insigni- 
ficant where  the  plasma  is  optically  thin. 

It  is  difficult  to  apply  this  method  to  plasmas  emitting  complex 
spectra  to  obtain  precise  values  of  n  .   In  addition  to  uncertainties 
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associated  with  inadequately  defined  atomic  constants,  there  are  problems 

in  defining  the  integrated-1 ine  intensity  I^     in  absolute  units  for  a 

^'total 
specific  transition.  Many  practical  considerations  in  the  measurement  of 

1      are  discussed  in  Ref.  36,  but  procedures  for  spectral -1 ine-wing 

total 
overlap  in  complex  spectra  are  especially  interesting. 

Uranium  spectra  show  no  isolated  line  structure;  therefore,  we  must 

find  the  peak  magnitude  and  FWHM  that  can  be  used  to  define  a  corres- 

49 
ponding  Voigt  profile,   a  theoretical  representation  of  the  actual  line 

profile  including  wings.  This  area  can  be  calculated  analytically  to 

yield  a  good  estimate  of  the  integrated  spectral -1 ine  intensity. 

B.        'Ppessuve-TeriTperature  Cox'relation    (PTC)    Technique 

Partial  pressure  estimates  of  plasma  constituents  can  be  based  a 

50  51 
pseudoanalytic  approach  such  as  the  PTC  technique.   '    This  method  uses 

a  temperature  profile  correlation  between  experimentally  determined  and 
calculated  temperature  profiles.  The  experimental  profile  is  typically 
established  by  the  Boltzmann  plot  method,  whereas  the  calculated  profiles 
are  computed  using  the  relative  (norm)  temperature  method.  Calculations 
of  radiation  specie  number  densities  and,  hence,  partial  pressures  are 
inherent  to  the  computed  profiles.  A  family  of  calculated  temperature 
profiles  is  generated  to  be  parametric  in  the  plasma  total  pressure.  Be- 
cause the  experimental  profile  is  an  independent  measurement,  intersection 
of  this  curve  with  that  of  the  calculated  profiles  implies  (with  the  aid 
of  the  Saha  equation)  a  plasma  number  density  and  partial  pressure. 

This  method  has  the  uncertainties  found  in  applying  the  relative 
(norm)  temperature  method  (Sec.  II. 4C)  as  well  as  experimental  inaccura- 
cies inherent  in  the  Boltzmann  plot  technique  (Sec.  II. 4A). 


Full -width  at  half  maximum  of  peak  intensity  value. 
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Uranium  plasma  number  densities,  as  determined  by  the  absolute-line 
method,  are  often  uncertain.  Pressure-temperature  correlation  removes 
some  of  this  uncertainty  and  furnishes  supporting  evidence  to  experi- 
mentally determined  densities. 


III.      URANIUM  PLASMA   EXPERIMENT  HARDWARE 
III- 1 .      Uraniwii  Plasma  Generation 

For  this  investigation  uranium  plasma  was  generated  by  a  direct- 

51 
current  uranium  arc  constructed  by  Randol   at  the  University  ot 

Florida.  Many  original  features  were  retained;  system  details  are  in 

Ref.  51.  There  have  been  some  important  changes  to  the  original  system 

which  will  be  described  in  this  report. 

/L Ui'anium  Plasma  Containment  Cell 

Figure  III-l  shows  the  uranium  plasma  containment  vessel.  The 
stainless  steel  vessel  is  designed  to  withstand  safely  cover-gas  pres- 
sure to  100  atmospheres.  It  can  also  be  operated  in  the  vacuum  mode 
down  to  at  least  300  torr.  Contact  is  made  between  the  tungsten  cathode 
and  uranium  anode  by  remote  movement  of  the  cathode  with  a  pneumatic 
electrode-drive  cylinder.  Both  electrodes  are  water-cooled.  The  gas 
distribution  head  can  give  directional  flow  to  the  incoming  cover-gas 
near  the  arc  electrodes.  The  viewport  windows  are  sealed  from  both 
sides  for  pressure  or  vacuum  operation.  The  gas  inlets  also  serve  as 
the  pumping  ports  when  an  evacuated  chamber  is  desired. 

IL iSig  qmsn  te  d  A  ssemb  ly 

Within  the  containment  vessel,  fastened  to  the  headplate,  are 
several  annular  water-cooled  copper  segments  (disks).  The  disks  are 
arranged  concentrically  around  the  anode-cathode  configuration  for  arc 
wall-stabilization.   Figure  III-2  shows  the  segmented  assembly  with  its 
orientation  to  the  electrode  configuration.  The  arc  column  length  is 
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shown  undersized,  and  the  plasma  column  was  actually  in  contact  with 
several  of  the  water-cooled  disks.  The  segmented  assembly  also  acts 
as  thermal  shield  for  the  pressure  seals  and  as  an  effective  particle 
shield  to  lessen  deposition  on  the  viewport  window. 

C .  Fvesino'e  Siistem 

A  schematic  representation  of  the  pressure  system,  designed  for 
flexibility  of  regulation  from  0-2000  psi,  is  shown  in  Fig.  III-3. 
Adequate  cell  exhaust  filtering  removes  uranium  which  might  escape  to 
the  atmosphere.  The  vent  on  the  downstream  side  of  the  cell  or  the 
roughing  pump  on  the  upstream  side  provides  particulate  venting.  The 
high-pressure  gas  supply  is  isolated  for  safety  in  case  of  electrical 
power  failure.  A  high-pressure  solenoid  valve,  normally  closed,  pro- 
hibits gas  flow  unless  the  solenoid  is  energized. 

D.  Coolant  System 

Enough  coolant  must  be  supplied  to  the  segmented  assembly,  anode, 
and  ballast  resistor.  The  cathode  is  primarily  convection-cooled  by 
the  surrounding  segmented  assembly.  The  60-psi  water-line  pressure  pro- 
vides adequate  cooling  for  an  arc  power  input  of  at  least  100  watts; 
however,  a  centrifugal  pump  provides  more  flow  if  necessary.  These 
features  are  shown  in  Fig.  III-4. 

/?._  Power'  Supply  System 

Two  650  A,  120  V  dc ,  diesel  motor  generators  arranged  in  a  series 
are  used  as  the  primary  source  of  electrical  power  for  arc  operation. 
They  can  be  operated  remotely  or  at  the  generator  controls  as  a  con- 
tinuously adjustable  voltage  supply.  Current  through  the  arc  circuit 
is  limited  by  air-  and  v/ater-cooled  ballast  resistors.  The  fuse  limit 
is  300  A,  250  V.  Current  is  adjusted  for  a  given  set  of  electrodes 
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during  arc  operation  by  parallel  switching  of  air-cooled  resistors  into 
or  out  of  the  circuit.  Typically,  this  operating  point  was  50  A,  20  V. 
The  uranium  arc  circuitry  is  shown  in  Fig.  III-5.  Current-voltage  mon- 
itoring is  accomplished  on  a  continuous  basis  by  Honeywell  stripchart 
recorders  (not  shown  on  Fig.  III-5). 
Ill- 2 . lira ni'um   Plai-nna  Stability 

Construction  of  a  dc  uranium  arc  was  not  difficult;  however,  de- 
velopment of  a  stable  dc  uranium  arc  which  would  allow  photoelectric 
diagnostics  demanded  extensive  effort.  For  spectroscopic  measurements, 
this  particular  uranium  plasma  was  required  to  be  very  stable  for 
density- temperature  measurements  and  marginally  stable  for  emission 
measurements.  Arc  current,  i  ,  was  used  as  a  stability  yardstick. 
Marginally  stable  implied  average  i  changes  <  -10%  over  the  entire 
length  of  data  collecting  time  with  a  maximum  of  ±  5%   for  instantane- 
ous i  changes.  Temperature-density  stability  is  indicated  by  similar 
average  1  limits  and  by  instantaneous  variation  of  1  5  2%.  The 
University  of  Florida  uranium  arc  evolved  from  a  free-burning  arc  to  a 
rather  sophisticated  wall-stabilized  arc  (Fig.  II1-6).  A  brief  account 
of  the  development  toward  increased  arc  stability  follows. 

/L FiMe-Bia-ninq  Ave 

For  simplicity,  a  free-burning  arc  under  a  static  helium  cover  gas 
was  used  (Fig.  III-6A).  Unfortunately,  motion  of  the  anode  and  cathode 
spots  was  Inherent  in  its  operation,  and  this  caused  unacceptable  move- 
ment of  the  arc  column  as  well  as  current-voltage  fluctuations.  There 
is  no  agreement  among  arc  physicists  as  to  the  reasons  for  these  spot 
movements;  however,  bibliographical  Information  can  be  obtained  from 
Ref.  46.  Many  techniques  were  applied  to  reduce  these  Instabilities, 
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such  as  sharply  pointed  and  polished  electrodes,  changes  in  arc  gap, 
changes  in  electrode  diameter,  and  different  cooling  rates,  but  all 
were  inadequate. 

'IL r,Vf;j-.9 tabilised  Ave 

Next,  a  gas  distribution  head  was  used  to  localize  large-scale 
niovenient  of  the  arc  column  (Fig.  III-6B),  and  cover  gas  (helium)  was 
directed  downward  from  the  gas  head  the  length  of  the  arc  column.  The 
flow  contact  with  the  anode  caused  a  divergence  of  the  cover  gas  at 
the  uranium  pellet  and  a  bell-shaped  cover-gas  flow  pattern  formed. 
The  cover-gas  flow  boundary  formed  the  "wall"  needed  for  arc-column 
local ization. 

C_. Tuhe-Stabilized  Avo 

Stabilization  was  enhanced  by  forced  containment  of  helium  gas 
flow  along  the  arc  column  (Fig.  III-6C).  A  quartz  tube  was  placed 
concentrically  with  the  electrode  vertical  axis  and  helium  was  forced 
the  length  of  the  quartz  tube.  The  restricted  flow  greatly  reduced 
the  helium  turbulence  and  its  effect  on  the  outer  and  inner  arc  column. 
However,  within  five  munutes  of  run  time  significant  particle  deposi- 
tion coated  the  quartz  tube  causing  unpredictable  intensity  attenuation 
and  prohibiting  long-duration  (photoelectric)  measurements. 

D. Segmented-Stabi lized  Ato 

The  segmented  assembly  (Figs.  III-2  and  TII-6D)  replaced  the  quartz 
tube  from  the  previous  case  and  reduced  the  deposit  problem  at  the 
viewport.  The  vacuum  (low-pressure)  segmented  arc  provided  stability 
for  photo-electric  intensity,  temperature,  and  density  measurements, 
while  the  high-pressure  segmented  arc  exhibited  marginal  stability 
acceptable  only  for  photo-electric  intensity  measurements.   The 
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segniented-arc  configuration  was  the  final  step  attempted  in  the  quest  for 
superior  arc  stability.  Even  with  marginal  stability  at  higher  pres- 
sures, this  system  was  employed  for  most  of  the  experiments  reported  here. 

III-?,. Pa  La  A  a,]  ids  Hi  on 

Spectroscopic  diagnostics  required  detection  and  analysis  of  radia- 
tion emitted  by  the  uranium  arc  plasma.  The  experimental  effort  of  this 
research  was  composed  of  two  broadly  defined  categories. 

(1)  the  measurement  of  intensity  (emission),  and 

(2)  diagnostics  for  temperature-density  determination. 

Figure  III-7  illustrates  necessary  components  for  simultaneous  measure- 

o 

ments  of  intensity  (to  2500  A),  temperature,  and  low-pressure  arc  den- 
sity. Intensity  measurements  which  extended  into  the  vacuum  uv  were 
performed  using  a  third  spectrograph  (McPherson  Model  218,  not  shown) 
designed  specifically  for  use  at  low  wavelengths.  Temperature  and  den- 
sity for  a  high-pressure  arc  were  inferred  from  photographic  spectral 

51        52 
analysis  completed  by  Randol   and  Mack   using  a  free-burmng  arc  at 

similar  I-V  conditions.  Details  are  in  Chapter  V. 

A.      Speabral  Line  Fvofile  cmd  Absorption  Data  for  the  Low-Pressure 
l\ra_ 

The  modified  BEM  (Sec.  II-4D)  was  used  to  determine  the  line-center 
absorption  coefficient  and  characteristic  plasma  temperature  for  the 
low-pressure  arc.  The  background  source  was  a  xenon  flashtube 
(EGG  FX-12-25).  The  firing  and  delay  schematic  is  in  Fig.  III-8.  Flash- 
tube  calibration  information  follows  in  Sec.  III-3C. 

Figure  III~7  shows  the  sequence  which  established  line  absorption, 
temperature,  and  density:  the  beam  splitter  passed  part  of  the  arc 
radiation  to  S-, .  A  rotating  refractor  plate  (quartz)  in  conjunction 
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with  spectrograph  (S-,)  used  as  monochrometer,  swept  a  particular  line 
of  interest  across  the  exit  slit  plane  of  S-, .  An  RCA  1P28  phototube 
placed  behind  the  S-,  exit  slit  responded  to  the  spectral  line  as  it 
moved  past.  That  line  profile  was  recorded  by  oscilloscopes  by  O-j  and 
0^.  At  the  instant  the  sweep  reached  the  central  wavelength  of  the 
line  profile,  the  xenon  flashtube  triggered.  This  enabled  the  flash- 
lamp  radiation  to  be  superimposed  on  a  specific  spectral  line  center 
(Fig.  II-3).  Flashtube  timing  was  accomplished  by  electronic  delay 
circuitry  shown  in  Fig.  III-8.  The  ground-state  particle  density  was 
determined  from  the  line  profile  of  oscillograms  O2  with  the  absolute 
line  method;  line  center  absorption  and  the  plasma  temperature  were  de- 
termined from  oscillograms  0-,. 

B. Fhotoeleati'-io  Intensity  Data 

o 

Photoelectric  intensity  measurements  were  made  from  2000  A  - 

o  00 

5500  A  for  the  low-pressure  uranium  arc  and  from  1050  A  -  6000  A  for 
the  helium-uranium  (high-pressure)  arc.  Intensity  of  the  low-pressure 
uranium  arc  as  a  function  of  wavelength  in  the  visible  and  near  uv  were 
recorded  by  using  part  of  the  data  acquisition  system  shown  in 
Fig.  III-2,,  A  scanning  spectrograph  $2  received  arc  radiation  reflected 
from  the  front  surface  of  a  beam-splitter.  The  phototube  response  was 
monitored  and  stored  digitally  by  a  signal-averager  [a  time-averaging 
digitizer  that  integrates  (smooths)  small  random  input  voltage  (arc) 
fluctuations]  which  resulted  in  ^^ery   reproducible  arc  intensity  traces 
as  a  function  of  wavelength.   Four  memory  areas  within  the  signal- 
averager  were  used  for  storage  of  the  spectral  intensity  I;^(;<),  where 
X  is  a  particular  LOS  position  in  an  arc  traverse.  A  four-point  Abel 
unfolding  for  spatial  resolution  of  the  arc  intensities  was  performed 
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to  obtain  the  spectral  emission  coefficient  e,(r),  where  r  is  the  radi- 
al  distance  from  the  arc  center.  The  limited  arc  stability  duration 
prevented  more  than  four  acceptable  LOS  intensity  measurements.  Oscil- 
loscope Op  was  used  to  trigger  the  signal-averager  sweep  at  the  desired 
wavelength.  Ultraviolet  transmission  through  the  optical  system  was 
carefully  investigated  to  be  certain  that  no  glass  was  present  and  to 

o 

study  the  system  attenuation  properties.  The  2537  A  Hg  spectral  line 
generated  by  a  very  stable  mercury  discharge  was  chosen  for  these  pur- 
poses. The  line  was  transiriitted  with  negligible  attenuation  and  con- 
tamination signals  produced  by  internal  reflections  within  the  spectro- 
graph Sp   were  removed  with  baffles. 

The  photomul tipl ier  used  with  S^  was  an  EMI-9514  with  a  sodium 
salicylate  window  which  acted  as  wavelength  shifter  from  the  ultraviolet 
to  the  visible.  The  phototube-sodium  salicylate  combination  greatly 
improved  the  system  wavelength  sensitivity  to  ultraviolet  and  vacuum 
ultraviolet  radiation.  Sodium  salicylate  was  ideal  for  use  in  the  uv 

and  vacuum  uv  because  it  possesses  a  nearly  constant  quantum  effi- 
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ciency  from  500  A  to  '-^   3300  A.    The  fluorescent  radiation  spectral 

o  o 

distribution  maximum  is  4300  A  and  ]0%   of  the  maximum  at  3800  A  and 

°   46  53  54 
5300  A,  '  '   which  conforms  to  the  maximum  wavelength  response  of 

many  photomul ti pliers. 

For  vacuum  uv  intensity  detection  the  McPherson  (Model  218)  spec- 
trograph designed  to  be  responsive  at  wavelengths  in  the  vacuum  uv  was 
used  in  conjunction  with  the  photomul tiplier  signal-averager  system 
previously  described.  This  particular  spectrograph  contained  magnesium 
fluoride-coated  (Al  +  MgF^)  optics,  a  2400-groove/mm  grating  biased  for 

0 

1500  A,  and  a  vacuum  capability  of  at  least  0.001  y.  The  vacuum  uv 
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region  is  reported  only  for  the  helium-uranium  arc  plasma  because  of 
its  relatively  higher  temperature  and  correspondingly  stronger  emission 
in  this  wavelength  region.  The  low-pressure  arc  did  not  have  noticeable 
emission  in  this  wavelength  region. 

In  some  instances,  it  was  physically  impossible  to  interchange  the 
source  with  a  reference  standard  at  the  same  location;  therefore,  we 
used  a  symmetric  arrangement  of  the  calibration  and  uranium  sources  in 
which  the  light-path  attenuation  for  either  source  was  identical.  We 
eliminated  unwanted  stray  radiation  from  within  the  spectrograph,  which 
is  particularly  critical  in  the  vacuum  uv  where  the  signal  detection  is 

o 

difficult.  Many  intensity-wavelength  scans  were  taken  below  2000  A 
with  a  deuterium  lamp  (Oriel,  C-42-72-12)  in  several  orientations  after 
each  helium-uranium  arc  run  to  insure  optimal  signal  transmission. 
These  signals  were  carefully  checked  to  minimize  higher  order  contami- 
nation and  internal  reflections. 

(\_  Ti itens i  bij  Ca  lihvation 

For  radiation  calibration  in  the  visible  and  near-uv  wavelength 
regions,  the  tungsten  lamp  and  the  positive  crater  of  a  carbon  arc  as 

standards  are  adequate;  however,  these  sources  are  unacceptable  for 

°  55 
lower  wavelengths  because  of  the  weak  intensities  below  2500  A.    In 

o 

fact,  below  2500  A  there  are  few  commercially  available  standard  cali- 
bration sources.  LTE  hydrogen  discharges  (fill  gas  is  either  hydrogen 
or  deuterium)  are  the  best  potential  sources,  but  they  require  exten- 
sive investigation  for  their  own  respective  properties.  The  theoreti- 
cal description  of  the  hydrogen  atom  is  essentially  complete,  and  once 
the  electron  densities  are  known,  the  intensities  are  computed  and 
cross-checked  by  experiment  to  provide  calibration  information  for 
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such  a  discharge.  Therefore,  two  calibration  sources  were  used:  a 
tungsten  lamp  for  the  visible  and  near  uv,  and  a  deuterium  discharge 
for  the  uv  and  vacuum  uv. 

Two  tungsten  filament  lamps  made  by  the  Eppley  Company  were  used 
for  intensity  calibration  in  the  visible  and  near  uv.  One  was  cali- 
brated by  the  NBS  and  designated  the  "standard  lamp";  the  second  served 
as  a  "reference  lamp."  The  reference  lamp  was  used  as  the  experiment 
standard  but  periodically  cross-calibrated  to  the  standard  lamp.  The 

calibration  curve  for  the  standard  lamp  is  given  in  Fig.  III-9.  The 

56 
associated  accuracy  of  the  values  was  stated  to  be  ±10%. 

A  deuterium  discharge  was  used  to  calibrate  intensities  below 

o 

2500  A.  The  lamp  had  a  suprasil  fused-silica  window  with  a  50%  trans- 

o  o 

mission  point  at  about  1750  A,  and  a  calibration  point  at  1662  A  was 
the  apparent  lower  wavelength  limit.  Figures  III-IO  and  III-ll  show 
the  wavelength  dependence  of  the  spectral  radiance  at  two  current  modes. 
The  absolute  intensities  as  calibrated  from  this  lamp  carry  a  ±10%  un- 
certainty,   verified  by  cross-check  of  the  deuterium  lamp  intensity 
against  the  tungsten  filament  standard  at  four  different  wavelengths 

o 

above  2500  A.  The  percentage  difference  between  the  quoted  deuterium 
intensities  and  the  cross-checked  intensities  was  always  within  the  un- 
certainty limits.  The  percentage  difference  tended  to  increase  toward 

o 

the  lower  vavelengths  and  the  calibration  points  below  2000  A  are  asso- 

o 

ciated  with  an  unknown  maximum  uncertainty  less  than  ±10%  up  to  1750  A. 

o 

The  last  deuterium  calibration  point  resides  at  1662  A  as  dictated  by 

o 

the  fused-silica  window  cutoff  of  the  discharge  lamp.  The  1662  A  cali- 
bration point  is  below  the  fused-silica  50%  transmission  wavelength  of 

o  57 

1750  A  and  is  quite  uncertain.    Calibration  for  intensity  data  to  the 
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LiF  cutoff  '(lOSO  A)  called  for  an  approach  which  did  not  require  expli- 
cit standard  calibration  points  at  the  lower  wavelengths. 

The  branching-ratio  technique  was  considered  for  absolute  intensity 
calibration  in  the  vacuum  uv  and  it  is  based  on  the  intensity  partition 
of  at  least  two  spontaneous  optically  thin  transitions  with  a  common 
upper  level. ^^'^^  The  intensity  ratios  of  these  two  lines  is  a  strict 
function  of  the  transition-probability  ratio  and  independent  of  plasma 
inhomogeneities  and  LTE.  To  apply  this  technique  in  the  vaccum  uv,  one 
transition  must  appear  in  the  visible  and  the  other  in  the  vacuum  uv. 
This  method  is  sound  provided  isolated  spectral  lines  exist  and  well- 
defined  transition  probabilities  are  used.  The  application  of  branching- 
ratio  calibration  to  uranium  spectra  fails  in  both  respects. 

Back-extrapolation  calibration  (BEC),  an  approximate  technique,  was 
used.  Helium-uranium  arc  intensity  data  were  collected  in  bandwidths 
900  A->1750  A  ->-  4300  A.  Intensities  above  1750  A  were  calibrated 
against  the  deuterium  discharge  and  tungsten  standard;  whereas  arc  in- 

o 

tensities  at  v^avelengths  below  1750  A  were  assigned  absolute  units  by 
BEC.  (Intensity  signals  from  the  low-pressure  uranium  arc  required 
calibration  down  to  2500  A  which  was  done  with  the  tungsten  standard.) 

Intensity  calibration  using  a  standard  source  is  essentially  a  sys- 
tem response  comparison  of  (in  this  case  a  spectrograph-photomultipl ier 
combination)  a  known  photon  flux  to  a  photon  flux  of  an  unknown  source. 
The  absolute  intensity  of  the  unknown  is  then  related  to  that  of  the 
standard  by: 


^ARc^'^)  -  ^std(^) -^7,7cn'  VgT^iy  -  ^"^  ^ 
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2     ° 
where  I„pp{A)  =  arc  intensity  [ergs/cm  -sec-A-str]  at  A, 

2    ° 
I   (a)  =  standard  source  intensity  [ergs/cm  -sec-A-str]  at  A, 

{J)   (A)  =  optical  path  reduction  factor  {%)   for  the  standard 
source  at  A, 

(j)„„p(A)  =  optical  path  reduction  factor  (%)  for  the  uncalibrated 
source  at  A, 

\l       [l]   =  response  signal  (volts)  of  optical  system  to  standard 
photon  flux  at  A,  and 

V,,np(A)  =  response  signal  (volts  of  optical  system  to  uncali- 
brated photon  flux  at  A. 

If  the  standard  and  uncalibrated  sources  are  in  identical  orientations, 
the  optical  path  reduction  ratio  (|-c,,pj(A)/cj)„np(^)  will  cancel.  The 
system  response  will  inherently  be  accounted  for  in  the  ^nori^)   ^"^d 
\l       (X)   signals.  Equation  (III-l)  is  valid  as  long  as  Ij-^„(A)  is  known, 

o 

in  this  case  to  the  1750  A  cutoff  of  the  deuterium  lamp  fused-silica 

o 

window.  To  apply  Eq.  (III-l)  to  wavelengths  less  than  1750  A,  a  system 

o 

response  below  1750  A  for  the  standard  source  intensity  as  a  function 
of  wavelength  must  be  assumed. 

o 

The  wavelength  dependence  of  system  response  below  1750  A  was 
assumed  flat  and  equal  to  1.0,  which  implies  that  the  combined  effect 
of  the  incident  (to  the  optical  system)  photon  flux  on  transmittance 
caused  by  losses  in  all  of  the  intercepted  optical  elements  is  negli- 
gible. This  is  categorically  not  the  case.  However,  such  a  tactic 
provides  a  straightforward  approach  to  a  conservative  estimate  of  the 
correct  absolute  intensities.  Relating  this  idea  to  Eq.  (III-l),  a 
unique  intensity  implies  a  unique  system  voltage  with  no  wavelength 
dependence.  Knowing  the  absolute  intensity  of  the  standard  and  its 
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o 

corresponding  voltage  at  1750  A  in  addition  to  the  arc  intensity  system 

o 

response  voltage  for  all  wavelengths  of  interest  below  1750  A,  the 
absolute  arc  intensity  may  be  estimated  by  solving: 


I   (1750  A)  V   (A) 
I   (A).-SID ARC _  (jjj_2) 

"^^  VgyplUSO  A) 


Vnnp(^)  retains  the  system  response  to  the  arc  intensity,  but  the 
wavelength  dependence  of  the  system  response  to  the  standard  intensity 
cannot  be  recorded  because  of  the  window  cutoff.  Therefore,  assuming 
the  negligible  optical  loss,  the  arc  intensities  are  underestimated  by 
the  composite  optical  loss  factors  and  whatever  molecular  absorption 
is  present. 

Whenever  an  approximation  is  applied,  as  was  the  case  with  the  in- 
tensity calibration  in  the  vacuum  uv,  an  indication  of  uncertainty  is 
valuable.  In  this  investigation  the  uncertainty  was  a  function  of  the 
attenuation  of  photon  flux  incident  on  the  optical  elements  of  the 
system  composed  of  a  10-mm  LiF  window,  two  spectrograph  mirrors,  the 
sodium-sal icylate  phototube  window,  and  the  grating.   (The  spectrograph 
mirrors  and  grating  were  Al-MgF2  coated).  A  composite  wavelength- 
dependent  transmission  curve  can  be  constructed  if  transmission/reflec- 
tance data  for  each  optical  component  are  available.  This  total  curve 
can  then  weight  the  appropriate  system  response  voltage,  thereby  re- 
covering (to  a  first  approximation)  the  photon  losses  as  the  field 
passes  through  the  optical  system. 

Reflectance-wavelength  information  is  available  in  the  litera- 
ture '   for  coated  and  uncoated  optics  as  is  transmission  data  for 
LiF.   '    The  grating  efficiencies  for  the  McPherson  2400-groove/mm 
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o 

grating  blazed  at  1500  A  (Al-MgF^  coated)  were  furnished  by  Quartz 
et  al .    Table  III-l  is  a  compilation  of  the  relevant  information  which 
allows  construction  of  a  composite  system  transmission  curve.  These 
data  provide  an  approximation  of  intensity  calibration  error  (using 
back-extrapolation)  based  on  system  optical  losses.  The  quantum  effi- 
ciency of  the  sodium  salicylate  phototube  window  was  assumed  to  be 

.  65   •  ,       ,   ,     n   J.I  ,    J     53,65,66 
approximately  unity   with  a  constant  wavelength  dependence. 

The  composite  efficiencies  of  Table  III-l  indicate  that  back- 
calibration  intensity  values  can  be  in  error  at  the  longer  wavelengths 

o 

(1750  A)  by  a  factor  of  6  and  more  toward  shorter  wavelengths.  A  cali- 

o 

brated  "signal"  detected  below  1200  A  is  suspect  because  of  the  sharp 
LiF  cutoff  near  1200  A  as  shown  by  Table  III-l. 

Use  of  the  xenon  flashtube  for  temperature-absorption  diagnostics 
required  a  calibration  for  dependence  of  brightness  temperature  on 
wavelength.  This  was  accomplished  by  comparison  (at  a  desired  wave- 
length) of  the  system  response  to  the  xenon  discharge  and  to  an  NBS- 
calibrated  tungsten-filament  lamp.    A  xenon  lamp  brightness  tempera- 
ture of  6745  K  ±  100  K  was  established  in  a  wavelength  range  of 
3600  A  ^  5400  A. 
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IV.      LOW-PRESSURE  URANIUM  ARC  DATA   REDUCTION 
IV-l.      Spectral   analjisis 

Spectr-a  emitted  by  the  low-pressure  uranium  arc  plasma  {'^^   300  torr) 
were  identified  manually  and  by  computer.    They  exhibited  a  mixture 
of  atomic  (UI)  and  singly  ionized  (UII)  uranium  line  structure,  but 
dominated  by  atomic  uranium  emission  lines.  For  diagnostic  purposes 

o  o  o 

the  3653.21  A  and  3659.16  A  UI  emission  lines  were  used.  The  3653.21  A 
line  is  reported   to  have  upper  and  lower  levels  at  31  166  cm   and 
3801  cm~  ,  respectively.   It  was  selected  primarily  because  of  its 
relatively  high-lying  lower  and  upper  levels,  which  enhanced  the  possi- 
bility of  evaluating  energy  states  populated  above  the  thermal  limit  if 

o 

partial  LTE  existed.  The  3659.16  A  line  was  chosen  for  the  same  diag- 

o 

nostic  purposes  (temperature  and  density  measurements)  as  the  3653.16  A 
line.  Its  transition  is  between  levels  27  941  cm"  and  620  cm  . 
The  first  excited  level  of  atomic  uranium  is  620  cm   and  it  would  be 
less  likely  to  meet  LTE  criteria.  A  comparison  of  temperature-density 
diagnostics  using  the  lines  with  two  distinctly  different  lower  levels 
served  as  one  useful  indicator  of  the  validity  of  LTE  assumptions. 
While  energy  level  was  a  prime  line-choice  factor,  the  spectral  lines 
were  also  selected  on  a  basis  of  wavelength,  relative  isolation  (from 
the  other  nearby  lines),  identification  certainty,  and  availability  of 
transition  probabilities  for  transitions  considered. 
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IV-2.      Temperature  Measwement 

The  electron  temperature  was  measured  by  the  modified  brightness- 
emissivity  method  (Sec.  II. 4).  First,  it  v/as  necessary  to  obtain  line- 
center  absorption  coefficients  as  a  function  of  arc  radius  by  applying 
flashtube  absorption  diagnostics  (Sec.  III.3A).  Line-center  absorption 

o  o 

profiles  using  the  3653.21  A  and  3659.16  A  UI  transitions  are  shown  in 
Fig.  IV-1  and  represent  absorption  measurements  taken  during  the  best 
conditions  of  arc  stability  and  flashtube  firing.  The  two  reported 
profiles  indicate  maximum  absorption  at  the  arc  center  which  is  usually 
the  region  of  maximum  temperature. 

Uncertainty  in  absorption  coefficient  values  was  both  numerical 
and  experimental.  The  unfolding  scheme  converting  k   (x)  ^  k   (r)   was 

^0        0 

responsible  for  the  numerical  uncertainty,  and  the  accuracy  of  the 
oscillogram  voltages  defined  the  primary  experimental  error  (k  indi- 
cates the  line-center  absorption  coefficient).  The  oscillograms 
contained  recorded  voltages  representing  the  flashlamp  intensity  atten- 
uated by  the  plasma  as  well  as  the  true  plasma  intensity.  The  plasma 
voltage  uncertainty  was  mainly  caused  by  arc  intensity  change  which  fell 
within  +2%   during  data  collection.   Flashlamp  attenuation  uncertainty 
resulted  from  arc  intensity  fluctuation,  flashlamp  intensity  variation, 
and  arc  cell  viewport  attenuation. 

Since  the  flashlamp  was  located  behind  the  arc  cell,  its  light 
output  passed  through  two  viewport  windows  in  the  segmented  assembly 
(Figs.  III-l  and  III-2);  whereas  the  arc  photon  field  passed  through 
only  one  on  its  path  to  the  detector.  During  arc  operation,  deposits 
occurred  on  both  viewports;  line-center  absorption  of  the  light  emitted  by 
the  flashlamp  was  recorded  photoelectrical  ly  at  four  points  from  the  arc 
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center  outward,  thus  the  values  at  further  distance  from  the  arc  center 

47 
were  more  uncertain.  Usher   has  considered  all  of  these  possible 

contributions  to  absorption  profile  uncertainty.  Using  an  argon  arc 

47 
with  the  same  experimental  apparatus  as  this  uranium  arc.  Usher   found 

a  fairly  constant  absorption  profile  uncertainty  across  the  arc.  This 
was  not  true  of  the  uranium  arc  because  the  viewport  deposition  prob- 
lems resulted  in  a  larger  absorption  coefficient  uncertainty  as  outer 
regions  of  the  arc  were  approached.  The  flashlamp  intensity  was  de- 
graded because  its  radiation  field  passed  through  two  viewpoints  (as 
opposed  to  one  for  the  arc  plasma  emission),  resulting  in  overestimates 
of  the  true  absorption  coefficient  values  at  all  points  along  the  arc 
radius.  Error  analysis  showed  that  these  estimated  uncertainties  were 
"^   10/0  at  the  centerline  to  "^  62°^  at  the  outer  point. 

The  apparent  uncertainties  of  radial  absorption  coefficients  in- 
fluenced the  temperature  computation  uncertainty.  Again,  the  error 

47 
analysis  of  Usher   was  used  to  obtain  temperature  uncertainties. 

Figure  IV-2  illustrates  the  radial  temperature  dependence  for  the  low- 
pressure  arc  plasma.  The  uncertainty  limits  generally  increase  as  the 
arc  center  is  approached.   If  we  consider  that  the  uncertainty  calcu- 
lated for  the  outer  radial  location  is  propagated  to  the  central 
locations,  it  is  possible  that  the  innermost  temperature  will  have  an 
associated  uncertainty  greater  than  that  of  the  outermost  temperature 
value.  Therefore,  while  the  outermost  temperature  had  the  largest 
experimental  error,  the  central  temperature  had  the  largest  total 
error.  Several  temperature  profiles  were  determined  using  both  tran- 
sitions (3653.21  A,  3659.19  A),  and  Fig.  I\/-2  shows  results  for  each 
transition.  The  3653  A  diagnosis  typically  showed  a  very  reproducible 
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o 

temperature  profile;  whereas  the  3659  A  transition  indicated  less  re- 
producibility. The  profiles  from  both  transitions  compared  rather  well 
and  offered  supporting  evidence  for  the  existence  of  LIE  conditions 
within  the  bulk  of  the  arc  column  down  to  fairly  low-lying  levels. 

Of  particular  interest  is  the  rather  flat  radial  temperature  de- 
pendence. There  have  been  numerous  flat  temperature  profiles  reported 

O  O  _  ')  Q 

for  various  wall-stabilized  arc  plasmas.      Plasma  confinement  using 
a  cooler  material  wall  generally  results  in  a  reduction  of  the  tempera- 
ture gradient  from  the  arc  center  toward  the  wall.  In  many  cases, 
conduction  loss  is  not  the  only  contributor  to  the  profile  curvature. 
Conductive  heat  transfer  along  the  arc  temperature  gradient  to  the  wall 
and  radiative  transfer  losses,  which  are  strongly  dependent  on  the 
plasma  conditions  and  type  of  radiating  species,  account  for  most  of 
the  energy  loss  from  an  arc.  That  a  flat  temperature  profile  in  wall- 
stabilized  arcs  can  be  realized  is  substantiated  in  the  literature  and 
not  in  controversy;  at  issue  is  the  rationalization  of  such  a  tempera- 
ture profile  for  a  wall-stabilized  uranium  arc  which  is  measured  by 
this  investigation. 

The  temperature  profile  for  wall-stabi 1 ized  arcs  must  be  consis- 
tent with  the  energy-balance  equation  given  by 


a  E^  =  V  •  F  +  V  •  Fn 
c       R 


(IV-1) 


where 


electrical  conductivity, 
electric  field, 
conduction  flux  density,  and 
radiation  flux  density. 
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The  implication  is  that  the  electrical  energy  input  to  a  steady-state 
arc  plasma  balances  the  radiation  and  conduction  losses  (convection 
being  neglected).  Conduction  loss  is  governed  by  the  temperature  pro- 
file gradient;  whereas,  radiative  losses  are  defined  by  the  equation 
of  radiative  transfer,  Eq.  (II-3). 

Temperature  profile  curvature  will  be  a  function  of  the  relative 
magnitude  of  the  conduction  or  radiation  loss  terms.  If  radiation 
effects  are  negligible,  conduction  becomes  the  dominant  loss  mechanism, 

resulting  in  an  approximately  parabolic  temperature  profile  with  a 

34 
comparatively  high  central  temperature  as  shown  in  Fig.  IV-3.    When 

radiation  losses  become  important,  the  central  temperature  is  generally 
lowered  and  the  profile  shape  flattened.  Also  included  in  Fig.  IV-3 
are  the  effects  of  self-absorption  on  the  temperature  profile.  In 
general,  having  zero  absorption  implies  a  smaller  central  temperature 
caused  by  the  larger  radiation  loss  term.  A  flatter  profile  also  re- 
sults because  most  of  the  curvature  is  caused  by  conduction  near  the 
wall.  When  absorption  is  included,  the  result  is  an  increase  in  cen- 
tral temperature  and  more  radial  curvature.  The  curvature  of  the 
profile  at  the  arc  center  is  typically  controlled  by  radiation  losses, 
while  at  the  arc  boundary  such  curvature  is  usually  controlled  by  con- 
duction losses. 

Several  parameters  affect  the  magnitude  of  conduction  and  radia- 
tion loss  terms.  Generally,  wall-stabilized  arcs  exhibit  high  radia- 
tion losses  at  higher  pressures,  greater  temperatures,  smaller  radii, 
and  at  greater  emission  density  of  the  dominant  radiating  plasma 
constituent.  "  '  '  '  The  current  and  pressure  of  the  low-pressure 
uranium  arc  would  indicate  that  the  temperature  profile  shape  would 
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not  be  highly  flattened  except  for  the  significant  radiation  flux  loss 
produced  by  such  a  dense  emitter  as  uranium.  The  striking  lack  of 
temperature  profile  curvature  (till  the  arc  boundary  is  approached)  can 
be  reasonably  attributed  to  the  dense  core  radiation  characteristic  of 
excited  uranium  in  a  wal 1 -stabil ized  configuration.    Shumaker   has 
shown  tliat  nitrogen  and  argon  wall -stabil  ized  arcs  operating  at  similar 
current,  voltage,  and  pressure  conditions  produced  temperature  profile 
shapes  of  different  degrees  of  curvature.  The  radiating  species  with 
the  denser  characteristic  radiation  (argon)  produced  the  flatter  tem- 
perature profile. 

From  LTE  considerations  the  normal  temperature  (Sec.  II. 4B)  for 
several  UI  transitions  was  calculated  to  be  about  5000  K  at  0.001  atm 
uranium  pressure.  The  5500  K  measured  temperature  indicated  a  mixture 
of  UI  and  UII  radiating  species  with  UI  dominating,  which  was  consistent 
with  the  spectral  analysis.   (See  Appendix  A  for  Saha  equation  number 
density  and  normal  temperature  curves).  The  probability  that  partial 
LTE  exists  also  implies  that  the  gas  temperature  differs  from  the 
electron  temperature--magnitude  of  difference  is  questionable.  This 
problem  has  been  studied  by  Gurevich   who  used  mercury  and  argon 
discharges  at  total  pressures  from  0.1  to  1  atm.  Their  technique 
essentially  monitors  the  distinct  cooling  of  the  electrons  to  the  gas 
temperature  and  subsequent  general  cooling  of  the  plasma.  The  intensi- 
ty output  tracks  the  cooling  and  can  be  followed  with  a  scope-phototube 
combination.  The  method  is  sensitive  to  electron-gas  temperature  dif- 
ferentials as  small  as  0.5^'.  The  conclusions  are  that  the  temperature 
differential  is  a  strong  function  of  total  (and  electron)  pressure  and 
weak  function  of  arc  current.  The  extrapolated  temperature  differential 
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for  the  present  arc  could  be  as  much  as  100  K  to  200  K  because  of  the 
partial  LIE  state  arc  plasma. 
IV- 3.      Density  Measiipement 

The  ground-state  particle  density  (UI)  for  the  low-pressure  urani- 
um arc  was  determined  by  the  absolute  line  method  (Sec.  II. 5A).  This 
method  involved  the  shape  definition  of  the  desired  spectral  line,  its 
absolute  intensity  calibration,  an  independent  temperature  measurement, 
and  uncertainty  estimation.  The  exact  line  profile  could  not  be  de- 
fined because  of  the  large  line-overlap  characteristic  of  uranium  spec- 

49 
tra.  A  Voigt  analysis   was  used  to  approximate  the  line  area  of  the 

o  C) 

same  lines  (3653  A,  3659  A)  used  for  the  temperature  measurements.  Be- 
cause these  data  were  taken  simultaneously  with  the  intensity- 
wavelength  information,  significant  temporal  fluctuation  should  be  com- 
mon to  all  and  treatment  of  fluctuation  was  unnecessary.  Line  intensity 
calibration  was  performed  with  a  tungsten-filament  NBS-cal ibrated  lamp. 
No  attempt  was  made  to  establish  a  radial  density  profile  because  asso- 
ciated errors  negated  the  effort.  Line-center  absorption  was  also 

accounted  for  by  appropriate  build-up  factors.  The  calibrated  UI  ground- 

14   -3 
state  density  for  the  low-pressure  uranium  arc  was  '^  7  x  10   cm  . 

This  value  is  an  order-of-magnitude  estimate  of  an  approximate  density 

radially  through  the  arc  column.  Assuming  uncertainties  in  spectral 

line  area,  gA,  and  temperatures  of  ±15,  ±50,  and  ±20%,  respectively, 

resulted  in  a  density  range  of  1.28  x  10   ->  7.89  x  10   cm"  .  This 

indicated,  from  Saha  analysis,  that  a  nominal  value  for  the  uranium 

total  pressure,  rounded  to  the  nearest  integral  logarithmic  pressure, 

would  be  -  0.001  atm  (1.3  x  10^^  cm"^). 


64 
The  electron  density,  n^,  for  the  low-pressure  uranium  arc  will 
generally  follow  the  curve  for  singly  ionized  uranium.  Saha  equation 
calculations  of  n^  for  a  plasma  with  a  total  uranium  pressure  of 

0.001  atm,  T  -  5500  K,  predicts  n^  ^  6  x  10^^  cm~^.  This  is  well 

11   -'^ 
within  the  LTE  ladder  criterion  of  3.8  x  10   cm   established  in 

Sec.  II. 1.  Measurement  of  the  electron  density  by  an  independent 
technique  such  as  a  line  broadening  analysis  was  too  uncertain,  pri- 
marily because  uranium  lines  lacked  isolation  and  definition.  There- 
fore, this  effort  was  abandoned  and  no  direct  experimental  evidence  of 

the  magnitude  of  the  electron  density  was  obtained.  However,  similar 

75 
work  of  Voigt   on  a  5500  K  uranium  arc  reported  n  to  be 

13   -3 
'^  5  X  10   cm  .  This  substantiated  the  assumption  of  n  large  enough 

for  the  present  arc  to  be  characterized  by  at  least  partial  LTE  nearly 

at  the  ground  level . 

One  of  the  more  subtle  aspects  of  this  density  determination  was 

the  implicit  assumption  of  complete  LTE.  Strictly  speaking,  this  is 

^ery   difficult  to  realize,  and  one  usually  resorts  to  a  commitment  of 

partial  LTE.  If  partial  LTE  is  valid,  an  exact  density  measurement 

(using  a  diagnostic  method  which  relies  upon  complete  LTE,  such  as  the 

absolute  line  method)  is  not  possible,  and  any  attempt  will  be  in  error 

because  of  this  apparent  conflict.  The  magnitude  of  this  error  can  be 

calculated  for  hydrogen  and  estimated  for  helium.  Unfortunately,  for 

the  case  of  uranium  only  a  qualitative  description  of  error  direction 
is  valid. ^^'77 

The  effect  of  the  complete  LTE  assumption  when  there  is  only  par- 
tial LTE  can  be  understood  by  considering  the  population  densities  of 
a  simple  energy  level  diagram  as  illustrated  in  Fig.  I\/-4. 
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Fig.  IV- 4.   Energy-level  diagram  with  thermal 
limit. 


When  complete  LTE  exists,  the  population  density  of  each  level  is  de- 
fined by  the  Boltzmann  Factor,  n  /n .  a  exp[-E  /kT].  However,  if 
partial  LTE  prevails,  the  population  of  levels  below  a  thermal  limit 
are  influenced  by  radiative  de-excitation  from  upper  levels  and  self- 
absorption,  particularly  with  transitions  terminating  at  the  ground 
level.  Detailed  definition  of  level  population  for  uranium  is 
virtually  impossible,  but  in  most  cases  one  can  argue  that  if  partial 
LTE  holds,  the  levels  below  the  thermal  limit  will  be  over-populated. 
Thus,  by  assuming  partial  LTE  and  using  complete  LTE  relations  to 
determine  the  ground-state  density,  an  underestimate  results,  probably 
proportional  to  some  function  of  the  thermal  limit  height  above  the 
ground  level.  The  thermal  limit  in  the  uranium  system  cannot  be  well 
established;  but  from  LTE  criteria  applied  to  the  low-pressure  uranium 
arc  (Sec.  II. 1),  it  is  likely  to  be  close  to  ground  level.   If  the 
thermal  limit  is  close  to  ground  level,  the  error  incurred  by  using 
complete  LTE  relations  to  determine  the  ground-state  population  is 
"small."  Because  this  is  an  order-of-niagni tude  measurement,  such  error 
is  most  likely  to  be  insignificant  at  these  temperatures  and  pressures. 
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TV- 4.      Emission  Coefficient  Determination 

Spectral  analysis  of  radiation  structure  emitted  by  the  lov;- 
pressure  uranium  arc  plasma  indicated  a  preponderance  of  known  UI  line 
spectra  mixed  with  some  identifiable  UII  structure.  The  species-mix 
ratio  is  not  available  because  of  uncertainties  involved  in  line 
identification  on  uranium  spectra.  However,  at  the  maximum  tempera- 
ture of  5500  K  a  Saha  ratio  of  "ijII'^'^UI  *^^"  approach  1.0,  and  cooler 
arc  regions  will  be  weighted  toward  larger  values  of  n^-r. 

Intensity  from  the  low-pressure  arc  was  observed  at  four  equally 
spaced  arc  positions.  Details  of  data  acquisition  are  in  Sec.  III.3B. 
The  spectrograph  used  for  this  study  was  a  Hilger-Ingis  (Model 
S05-10000)  modified  to  a  rapid-scan  capability  (1.28  s  from  2000  to 

O 

6000  A).  A  1200-grooves/mm  grating  was  used  with  a  blaze  angle  of 

o 

17.2°  corresponding  to  a  wavelength  of  5490  A  in  the  first  order.   In- 
tensity calibration  was  performed  using  a  tungsten-filament  NBS- 
calibrated  standard. 

Preliminary  experiments  indicated  a  rapid  decrease  of  arc  intensi- 

o  o 

ty  below  3500  A  and  very  little  signal  at  wavelengths  less  than  2500  A. 
Steinhaus  et  al .   '   also  indicates  little  structure  between  2000  and 

o 

2500  A  for  UI.  We  decided  that  emission  from  this  particular  arc 

o 

plasma  below  2000  A  would  be  relatively  insignificant,  and  the  effort 
required  to  detect  potential  vacuum  uv  emission  was  not  justified  for 
this  plasma. 

The  calibrated  arc  intensity-wavelength  data  is  plotted  in 
Fig.  I\/-5.   Four  similar  sets  of  intensity  data  were  generated  for  each 
arc  burn-,  the  data  plotted  in  Fig.  l\l-5   corresponds  to  the  central  arc 
region  only.  The  data  acquisition  method  (signal  averaging)  produced 
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very  reproducible  intensity-wavelength  data  sets;  only  major  arc  inten- 
sity fluctuations  could  not  be  removed.  When  these  major  fluctuations 
occurred  the  data  were  discarded,  thus,  the  data  plotted  in  Fig.  IV-5 

is  in  error  by  the  intensity  uncertainty  of  the  standard  tungsten  lamp 

55 
calibration  source,  which  is  approximately  ±10%.    The  calibration  of 

uranium  plasma  intensity  was  done  using  the  "Plasma  Source  Calibration 

78 
Program."  '     However,  intensity  data  sets  taken  in  the  outer  arc  re- 
gions are  in  error  by  as  much  as  +12  to  20%  because  deposition  on  the 
viewport  interfered  with  the  optical  path  from  the  plasma  source  to  the 
spectrograph  entrance  slit. 

The  four  intensity-wavelength  data  sets  were  spatially  resolved  by 

o 

a  four-point  numerical  Abel  unfolding  after  averaging  over  100  A  band- 
widths.  The  results  for  this  emission  coefficient  calculation  in  the 
central  arc  location  are  shown  in  Fig.  IV-6.  Error  sources  in  these 
data  are  attributed  to  the  absolute  intensity  calibration,  viewport 
deposition,  and  four-point  spatial  resolution.  The  calibration  and 
deposition  errors  have  already  been  mentioned  and  are  carried  over  to 
the  emission  coefficient  determination.  Spatial  resolution  error  is 
much  more  difficult  to  address.   It  is  composed  of  error  caused  by 
propagation  of  experimental  uncertainties  and  the  high  probability  of 
using  a  non-optimal  zone  number  for  the  unfolding  (statistical  error). 
Experimental  error  propagation  through  the  inversion  process  was 
estimated  by  unfolding  the  data  with  and  without  such  error  included. 

The  statistical  error  was  estimated  by  relying  on  the  analysis  of  Kock 

44 
and  Richter,   based  on  the  form  of  the  intensity  profile,  number  of 

rings  chosen,  smoothness  of  profile,  and  maximum  outer-zone  experi- 
mental error.  An  exact  match  did  not  exist  and  extrapolation  was 
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applied.  Following  the  above  approach,  the  error  for  the  central  zone 
emission  coefficient  shov;n  in  Fig.  l\l-6   should  fall  within  -14%  to  +23%. 
Credibility  of  the  shape  and  location  of  the  important  features  on 
the  wavelength  scale  is  questionable;  unfortunately,  there  are  no  other 

data  recorded  at  the  same  conditions  for  comparison.  However,  compari- 

18  19  79 
son  data  exist  '   '   for  qualitative  assessment  of  credibility. 

79 
Krascella   tried  to  approximate  the  intensity-wavelength  distri- 
bution of  an  Argon-UF,  system  at  temperatures  varying  from  5000  to 
9000  K.  Level  populations  were  calculated  through  equilibrium  rela- 
tions, whereas  partition  functions,  atomic  constants,  and  observed  line 

69 
location  were  extracted  primarily  from  Corliss  and  Bozman.    Using 

this  information  the  possible  integrated  line  intensities  were  computed 

o 

and  averaged  over  100-A  bandwidths.  This  semiempirical  technique  lacks 
information  about  important  quantities  such  as  statistical  weights, 
accurate  transition  probabilities,  uncertain  UF^  decomposition  schemes, 
and  offers  very  incomplete  uranium  line  structure  tables.  However, 
this  work  was  valuable  in  establishing  observable  spectral  distribution 
of  electromagnetic  radiation  from  uranium  plasmas  at  various  conditions 
and  was  useful  for  comparison  to  our  results  even  though  Krascella 
estimates  had  gaps.  The  comparison  at  least  substantiated  a  trend  in 
the  emission  coefficient  shape.  No  comparison  was  made  between  abso- 
lute values  because  of  the  diversity  in  systems. 

A  comparison  is  shown  in  Fig.  IV-7  where  the  Krascella  data  were 
shifted  in  magnitude  to  be  superimposed  upon  our  results  which  follow 
the  shape  trend  shown  by  Krascella.  The  irregularities  in  the 
Krascella  data  (caused  by  the  lack  of  experimental  line  data  in  the  gap 
regions)  are  expected  and  do  not  detract  from  the  conclusion  that  the 
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two  data  sets  are  supportive.  Thus,  we  have  established  corroborative 

evidence  of  the  shape  credibility  of  the  low-pressure  arc  emission  co- 
efficient. 

The  work  of  Steinhaus  et  al .   '  also  lends  support  to  the  major 

o 

peak  location  around  4000  A   (for  UI)   and   to  the  validity  of  the  rapid 

o 

emission  decrease  below  3500  A.  The  range  of  experimentally  observed 
UI  levels  extends  to  40  000  cm   on  the  energy-level  diagram  ('^  5  eV). 
This  implies  there  are  no  known  UI  transitions  below  approximately 

o 

2500  A.  Also,  the  bulk  of  known  lines  established  by  Steinhaus  for  UI 

o 

falls  within  3300  to  6000  A,  which  also  supports  our  information. 

The  Planck  function  for  a  blackbody  temperature  of  y   5500  K  is 
also  plotted  in  Fig.  IV-7  along  with  the  emission  data  (different 
units).  Since  the  magnitudes  of  the  arc  emission  and  Planck  intensity 
are  plotted  logarithmically,  the  difference  of  these  values  at  a  spe- 
cific wavelength  is  an  indication  of  the  absorption  coefficient,  k,. 

o 

The  value  of  k-,  is  actually  an  average  (over  100  A)  because  the  emis- 

A 

o 

sion  coefficient  is  also  averaged  over  100  A.     At  all  wavelengths 


^A 


0.034,   thus   implying  an  optically  thin  plasma.     Line-center  ab- 


sorption coefficients  can  be  significantly  greater  for  some  lines  as 

o  o 

illustrated  by  the  line-center  k.    for  the  3653  A  and  3659  A  UI  lines. 


V.         Ul'lLlim-UiiMlUM  ARC   'DATA   REWCTION 

The  objective  of  this  phase  of  research  was  the  investigation  of 
uranium  plasma  emission  at  higher  temperature  (than  that  of  the  low- 
pressure  uranium  arc)  from  the  visible  through  portions  of  the  vacuum 
ultraviolet  wavelength  region.  The  experimentally  observed  levels  of 
UII  (singly  ionized  uranium)  extend  up  to  50  000  cm  ,  ^''       and  the 
ionization  energy  is  about  100  000  cm  .  Hence,  UII  line  structure  is 
likely  to  appear  in  the  visible,  uv,  and  vacuum  uv  wavelength  regions. 

For  this  reason  a  uranium  arc  which  produced  strong  UII  radiation  and 

51 
operated  under  conditions  similar  to  those  observed  by  Randol   and 

Mack   was  used. 

The  uranium  plasma  for  the  present  study  v;as  generated  with  the 

identical  configuration  as  the  low-pressure  arc  (Figs.  III-l  to  5). 

Helium  cover-gas  at  3  atm  was  added;  the  arc  current  and  voltage  were 

maintained  at  30  A  and  35V,  respectively.  Current-voltage  and  arc 

emission  characteristics  were  controlled  to  duplicate  those  determined 

51 
by  Randol   to  infer  temperature  and  density  from  his  photographic 

diagnostics.  This  was  an  important  consideration  because  our  uranium 

arc  was  too  unstable  for  accurate  photoelectric  temperature  and  density 

diagnostics  such  as  those  performed  on  the  low-pressure  arc.  This 

inference  seemed  reasonable,  if  spectral  similarity  (line  location, 

half-width,  and  peak  values)  for  arc  plasmas  generated  by  cascade  and 

free  burning  systems  at  the  same  pressure-current-voltage  conditions 

could  be  established  within  the  error  limits  of  Randol 's  temperature 

and  density  measurements. 

73 


V-1.      Spectral  Analysis  74 

The  line  spectrum  emitted  from  the  helium-uranium  arc  was  identi- 
fied manually  and  cross-checked  with  the  line  identification  computer 
analysis  of  Kylstra.    Identified  lines  were  highly  correlated  with 
the  known  UII  spectrum.  Only  a  few  helium  line  possibilities  were  ob- 
served, which  indicated  a  relatively  pure  uranium  plasma.  This  spectrum, 
as  a  whole,  was  significantly  different  in  comparison  with  the  low- 
pressure  arc  plasma,  and  provided  some  evidence  that  operation  from  a 
low-  to  a  high-pressure  arc  was  accompanied  by  a  corresponding  shift  of 
radiation  dominance  to  UII  in  addition  to  a  shift  from  lower  to  higher 

temperature. 

51 
Comparison  of  helium-uranium  arc  spectra  sets  (ours  and  Randol's  ) 

was  first  done  on  a  wavelength  basis.  When  superimposed,  the  two  sets 
of  spectra  were  virtually  indistinguishable.  Both  line  sets  were  com- 
posed of  line  spectra  at  the  same  wavelength  locations  and  emitted  by 
the  same  plasma  constituent;  namely,  singly  ionized  uranium.  This  was 
necessary  but  not  enough  to  justify  the  assumption  of  similar  tempera- 
ture and  pressure  (T-P||tt).  The  other  necessary  factor  (provided  most 
of  the  lines  were  sharp)  to  insure  similar  T-P,,^j  conditions  was 
consistency  between  the  two  sets  of  peak  and  half-width  values.  De- 
tailed study  showed  little  discrepancy  for  the  lines  examined  at  the 

51 
arc  centerline  used  in  Randol's   diagnostic  analysis.  Further,  peak 

and  half-width  values  of  these  lines  compared  with  those  recorded  at 
different  current,  voltage,  and  pressure  conditions  showed  discernible 
differences  which  indicates  the  sensitivity  of  the  comparison.  There- 
fore, the  centerline  temperature  and  pressure  of  the  constricted 

51 
helium-uranium  arc  plasma  was  approximately  characterized  by  Randol 


75 

to  be  between  8000  and  9000  K  and  0.01  atni  UII  partial  pressure  rounded 
to  the  nearest  integral  logarithmic  pressure.  Because  of  the  inferred 

temperature  and  |:)ressure,  only  a  brief  overview  will  be  given  (Sec.  \/.2 

51 
and  v. 3)  of  tlie  methodology  used  by  Randol   to  determine  these  param- 
eters. 

F-g. Tempevatur'S  Msas urement 

The  central  temperature  of  the  helium-uranium  arc  plasma  operated 
at  30  A,  36  V,  and  3-atm  total  pressure  was  estimated  by  the  Fowler-Milne 
method  and  measured  by  the  relative  Boltzmann  plot  method.  The  normal 
temperature  for  this  plasma  was  estimated  to  range  between  7000  and 
9000  K  for  a  corresponding  pressure  range  of  0.001  to  0.1  atm  total 

o 

uranium  pressure,  using  the  4171  A  UII  transition.  Slight  off-axis 
maximums  in  the  emission  profile  indicated  the  probability  of  close 
proximity  to  the  normal  temperature.  To  define  temperature,  the  corres- 
ponding pressure  must  be  known  because  temperature  and  pressure  must  be 
consistent.  Because  the  spectrum  taken  from  this  plasma  showed  a  clear 
majority  of  UII  structure,  and  Saha  analysis  indicated  (at  0.01  atm)  a 
dominance  of  UII  particles  in  the  7000  to  9000  K  temperature  range,  the 
probability  of  close  proximity  to  the  normal  temperature  was  increased. 

The  centerline  temperature  was  then  measured  for  the  above  conditions 

51 
using  the  Boltzmann  plot  method  and  found  to  be  8113  K  t  8%.    The 

centerline  temperature  of  our  helium-uranium  arc  plasma  operated  at 
similar  conditions  was  therefore  inferred  to  be  ^   8000  K. 
V-Z.      PerisitLj  Measurement 

The  central  UII  partial  pressure  for  the  helium-uranium  arc  was 
measured  by  the  absolute  line  method  and  cross-checked  by  the  pressure- 
temperature  correlation  (PTC)  technique  (Sec.  II. SB).  This  diagnostic  was 
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similar  to  that  applied  to  the  low-pressure  arc.  Order-of-magnitude  dis- 
agreement occurred  between  the  absolute  line  values  and  those  determined 

51 
by  vaporization  and  PTC  studies.    In  fact,  the  absolute  line  values 

were  always  lower  than  the  others.  These  discrepancies  were  explained 
(to  some  degree)  by  assuming  an  inhomogeneous  pressure  profile  radially 
across  the  arc  and,  as  well,  the  characteristic  uncertainties  asso- 
ciated with  most  density  measurements.   In  any  event,  the  centerline  UII 
density  reported   for  a  3~atm  helium-uranium  arc  vjas  -  3.8  x  10   cm 
(0.042  atm)  at  a  central  temperature  of  'i-  8100  K.  By  Saha  analysis  this 
pressure  and  temperature  imply  a  total  uranium  pressure  of  '^  0.055  atm 
and  an  electron  density  of  4.3  x  10   cm  . 

The  temperature  measurement  assumed  partial  LTE;  however,  the  spec- 
troscopic density  determination  required  complete  LTE.  As  described  in 
Sec.  IV. 3,  assuming  complete  LTE  where  only  partial  LTE  was  assured 
caused  an  underestimation  of  the  ground-state  densities;  the  magnitude 
of  error  vias   a  function  of  the  thermal  limit.   For  the  helium-uranium 
arc  a  rather  high  electron  density  was  achieved,  which  easily  satisfied 
the  ladder,  and  possibly  the  hydrogenic  LTE  criterion  as  well.   Based  on 
these  considerations,  the  ground-state  underpopulation  factor  (if  it 
could  be  calculated)  would  be  small  compared  to  the  other  uncertainties 
found  in  the  density  measurement.  The  collective  error  applied  to  this 
density  measurement  implied  slightly  better  than  an  order-of-magnitude 
estimate. 
V-d.      Emission   Coeffiaient  Determination 

Radiation  from  a  helium-uranium  arc  operated  at  a  total  pressure  of 
3  atm  and  -^v  1000  W  arc  power  was  detected  without  any  collimating  lenses. 
The  arc  chamber  was  continually  purged  with  helium  and  the  vacuum 


11 

spectrograph  (McPherson  Model  218)  v;as  maintained  at  •\'   1  n.  Separating 
the  pressure-vacuum  (arc  cell-spectrograph)  interface  was  a  3/8-in. 
Lif  (lithium  fluoride)  window.  The  spectrograph  grating  contained 

o 

2400  grooves/mm  and  was  blazed  for  1500  A  in  the  first  order.  The 
sodium  sal icylate-phototube  combination  (described  in  Sec.  III.3B)  was 
used  to  detect  photons  transmitted  by  the  spectrograph. 

After  initial  detection  of  a  vacuum  uv  signal,  we  tried  to  elimi- 
nate possible  undesirable  signals  caused  by  internal  reflections,  stray 
light,  etc.  The  McPherson  spectrograph  had  a  history  of  reflection 
problems,  so  optical  blockouts  (baffles)  were  used  on  the  two  spectro- 
graph mirrors  and  the  grating.  The  geometry  of  the  baffles  allowed 
only  the  central  image  of  a  light  source  to  pass  through  the  spectrograph 
to  the  exit  slit,  and  the  spectrographic  internal  reflections  were 
el iminated. 

Arc  intensity  data  v/ere  collected  in  two  wavelength  segments,  from 
1050  A  to  1750  A  and  from  1750  A  to  4300  A.  These  bandwidths  were  used 
because  the  deuterium  calibration  standard  had  no  valid  calibration 

o  o 

values  lower  than  1750  A.  Intensities  at  wavelengths  above  1750  A  were 
calibrated  as  described  by  Eq.  (III.l).  The  calibration  standard  used 

o 

at  wavelength  >  1750  A  was  the  Oriel  deuterium  discharge.  This  region 
was  cross-checked  at  several  wavelengths  using  the  tungsten  filament 
standard  as  a  reference,  and  most  of  the  uncertainty  associated  with  the 
deuterium  calibration  values  was  removed. 

Accurate  calibration  of  intensity  data  found  in  the  vacuum  uv  region 

o 

below  1750  A  by  ordinary  methods  was  impossible.  Many  of  the  difficul- 
ties encountered  are  discussed  in  Sec.  III.3C;  application  of  back- 
extrapolation  calibration  (BEG)  was  necessary  (see  Sec.  III.3C). 


78 

Because  the  intensity  data  were  collected  without  a  collimating  lens, 
no  detailed  spatial  resolution  was  performed  on  the  observed  intensi- 
ties. However,  a  homogeneous  plasma  of  1 -cm  depth  (the  actual  arc 
plasma  diameter)  was  assumed.  The  intensity  was  converted  to  an 
approximation  of  the  arc-centerl ine  emission  coefficient  by  simply 
weighting  each  intensity  value  by  the  inverse  of  the  arc  plasma  depth 
(in  this  case  1/1).   In  general,  this  procedure  yielded  a  conservative 
estimate  of  emission  coefficient  and  a  slight  distortion  of  the  true 
shape. 

The  calibrated  and  converted  intensity  values  are  shown  in  the 
approximate  form  of  emission  coefficient  values  in  Fig.  V-1 .  There  is 

o 

moderate  emission  in  the  1000  to  2000  A  bandwidth  mostly  caused  by 
overlapping  line  structure.  This  bandwidth  has  been  rescaled  in  Fig.  V-2 
to  expose  more  spectral  detail.  Figure  V-2  indicates  notable  line 

o 

structure  and  one  distinct  emission  peak  at  '\^  1550  A,  partly  caused  by 
grating  response  at  the  nominal  blaze  angle. 

o 

Grating  scan  uncertainty  of  +10  A  and  rather  large  slit  widths  made 

precise  wavelength  location  of  the  line  structure  unobtainable  for  corn- 
on 
parison.  Kelly,   however,  lists  a  rather  intense  cluster  of  lines  at 

o         o  o 

1575  A,  1579  A,  and  1585  A--consistent  with  our  results.  The  emission 
data  shown  in  Fig.  V-2  are  one  of  many  sets  taken  (at  similar  arc  condi- 
tions) which  were  cross-referenced  to  eliminate  noise  and  insure  that 
most  of  the  residue  be  the  desired  signal. 

Figure  V-1  shows  the  remainder  of  the  spectrum  to  4300  A.  Many  of 
the  known  UII  and  some  of  the  stronger  UI  lines  vyere  potentially  iden- 
tified. Many  unidentifiable  lines  could  be  of  UII  origin.  There  are 
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o  o 

distinct  peaks  at  approximately  2300  A  and  2900  A.  The  theoretical 

o 

predictability  of  these  peaks  and  the  one  at  1500  A  will  be  addressed 
in  Chapter  VI. 

Vacuum  uv  signal  authenticity  and  order  contamination  caused  by 
overlapping  orders  were  investigated.  Light  emitted  by  the  uranium  arc 
was  passed  through  the  vacuum  spectrograph  system  and  detected  as  usual 
with  one  exception--a  test  material  composed  of  either  Lit",  quartz,  or 
glass  was  placed  just  behind  the  entrance  slit.  Arc  intensity  data 
collected  indicated  that  the  system  was  detecting  a  true  vacuum  uv 
signal,  and  that  contamination  from  other  orders  of  the  grating  was 
negl igible. 

Uncertainty  in  the  absolute  value  of  the  emission  coefficient  is 
a  function  of  several  sources,  such  as  calibration  of  the  standard 
radiation  source,  minimal  spatial  resolution,  standard  source  position- 
ing, digital  processing,  grating  scan  error,  and  arc  fluctuation.  For 

o 

higher  wavelength  (>1750  A)  the  emission  coefficient  value  error  is 
probably  within  j30%.  However,  it  is  much  more  difficult  to  identify 

o 

value  error  to  wavelengths  less  than  1750  A.  Back-extrapolation 
calibration  in  this  region  assumed  negligible  system  (optical)  losses 
which  is  not  the  case.  A  good  indicator  of  value  uncertainty  at  these 
wavelengths  is  the  system  efficiency  data  in  Table  III-l.  The  dashes 
in  Fig.  V-1  represent  the  uncertainty  limits  on  data  at  selected  wave- 
lengths in  the  vacuum  uv.  Clearly,  as  lower  wavelengths  are  approached, 
the  losses  begin  to  dominate  and  are  really  indicative  that  all  photons 

o 

emitted  below  1200  A  are  unlikely  to  be  detected.  The  error  in  emis- 

o 

sion  coefficient  values  at  wavelength  >  1750  A  does  not  have  this  loss 
component  because  the  calibration  is  more  precise. 
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o 

An  extension  of  these  data  to  6000  A  was  accomplished  by  comparison 
to  uranium  plasma  emission  coefficient  data   taken  using  this  arc  sys- 
tem and  is  shown  in  Fig.  V-3.  The  line  structure  in  this  figure  refers 
to  the  latest  measurements  at  3-atm  total  pressure;  the  dashed-line 
curves  represent  previous  results.  Tlie  deviation  in  magnitude  can  be 
explained  by  our  more  sophisticated  and  accurate  data-acquisition  system. 
Even  so,  the  comparative  magnitudes  remain  within  reason.  Comparative 
shapes  (along  the  wavelength  scale)  offer  a  high  degree  of  correlation, 
as  expected  if  past  and  present  measurements  were  made  properly.  These 

o 

factors  validate  the  extrapolation  of  the  present  data  to  6000  A  with 
substantial  credibility.  Thus,  a  spectral  emission  coefficient  is  now 
defined  for  a  uranium  plasma  (8000  K,  0.01-atm  uranium  pressure,  and 
1-cm  plasma  depth)  from  1200  to  6000  A. 

Figure  V-4  shows  the  present  results  with  the  corresponding  Planck 

function  and  other  comparable  measurements  reported  in  the  litera- 

77   ?3  ?7  83 

ture,  '  '   along  with  theoretical  predictions  made  by  Parks  et  al. 

For  the  moment  Parks'  results  will  be  accepted  and  their  validity  ex- 
amined in  Chapter  VI.  Because  the  graph  is  semilogarithmic,  Kirchoff's 
Law  provides  a  ready  means  to  estimate  directly  the  absorption  coeffi- 
cient by  merely  subtracting  the  value  difference  between  plasma  emission 

and  the  Planck  function  at  a  given  wavelength. 

22 
For  comparison  it  is  necessary  to  remember  that  the  Miller   and 

7 1 
Marteney  et  al .   experiments  used  UFg  as  the  discharge  gas;  the  Florida 

experiment  vaporized  metallic  uranium.  Figure  V-4  clearly  indicates  the 
differences  in  emission  coefficient  wavelength  dependencies  among  experi- 
ments. At  similar  plasma  conditions  and  compositions,  differences  in 
shape  should  be  minimal.  A  distinct  fall -off  in  emission  coefficient  is 
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apparent  between  2500  and  4000  A  in  both  the  Miller  and  Marteney  results 
that  is  not  evident  in  the  Florida  values.  This  discrepancy  may  be 
caused  by  strong  UF^  absorption  in  this  wavelength  region.  This  can  be 

supported  by  Fig.  V-B  which  is  a  reproduction  of  some  of  the  latest  UF^ 

84  85  8fi  87  88 

photo  absorption  '   and  electron  impact  cross-section  data     ' 

indicating  significant  UF^  absorption  at  these  wavelengths.   It  is  un- 
likely that  other  candidates  such  as  UFr,  UF^,  and  F  play  a  major  role 

in  light  absorption  within  UFg  discharges  because  their  mean  free  paths 

89  90 
are  typically  on  the  order  of  centimeters. 

Recent  and  fairly  conclusive  evidence  shows  that  Miller's  meas- 
urements may  exhibit  a  substantial  absorption  effect  from  cold  layers 
of  UFg,  particularly  in  the  2000-  to  4000-A  bandwidth.  Spector^^  per-' 
formed  a  low-temperature  (700  to  1400  K)  ballistic  piston  UFg  absorption 
experiment  with  results  remarkably  similar  in  magnitude  and  shape  to 
those  of  Miller  whose  temperature  was  a  reported  10  000  K.  The  implica- 
tion of  this  similarity  is  that  while  the  uranium  plasma  in  the  Miller 
experiment  may  have  been  at  10  000  K,  its  emission/absorption  charac- 
teristics were  masked  by  such  properties  of  UF^.        ■  ' 

Unfortunately,  the  Miller  data  do  not  extend  below  2500  A;  however, 
the  Marteney  and  Florida  results  do  extend  into  the  vacuum  uv.  Both 
indicate  a  small  emission  peak  between  1400  and  1800  A.  The  Florida 

o  o 

peak  appears  at  1500  A,  while  the  Marteney  peak  is  located  at  1650  A. 
The  wavelength  displacement  between  the  two  peaks  cannot  be  explained  by 
experimental  error  and,  therefore,  is  attributed  to  the  nature  of  each 
plasma  and  its  associated  emission/absorption  mechanisms.  The  small 

o 

peak  shown  by  the  Marteney  curve  at  'x^  1650  A  is  inconsistent  with  the 

o 

relatively  large  UFg  absorption  cross  section  in  the  2000  to  4000  A 
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bandwidth.  However,  work  by  Srivastava  et  al .  ,   McDiarmid,   and 


Trajmar   has  substantiated  significant  variation  in  the  UF^  absorption 

o 

cross  section  in  the  1500  to  2000  A  bandwidth  as  illustrated  in 

Fig.  V-S.      (For  absolute  units  the  common  point  between  Figs.  V-5  and 

o 

V-6  is  at  the  peak  value  near  2255  A.)  In  fact,  there  appears  to  be  an 

o 

absorption  window  at  1650  A  which  may  account  for  the  emission  peak  at 
that  wavelength  from  the  Marteney  measurement.  This  may  also  explain 
the  vacuum  uv  peak  shift  between  the  Florida  and  Marteney  results. 

Now  we  examine  UF^  photoabsorption  effects  on  the  present  arc  data 
and  the  Marteney  data  because  it  may  explain  some  of  the  discrepancies 
in  shape.  The  Florida  arc-emission  coefficient  data  shown  in  Fig.  V-1 
were  folded  with  several  UF^  photoabsorption  strengths  defined  by  layer 
thickness  (as  in  the  Marteney  experiment)  and  molecular  density.  The 
concentrations  correspond  to  UF^  to  approximately  0.01  to  0.1  atm  at 
room  temperature.   (The  details  of  the  folding  are  found  in  Appendix  E.) 
The  results  are  shown  in  Fig.  E-1  through  E-4.  They  indicate  a  rather 

o 

severe  emission  reduction  in  the  2100  to  2900  A  bandwidth,  a  possible 

o  o 

peak  emerging  at  1750  A,  and  the  original  peak  at  1500  A  reduced  and 
shifted  toward  the  lower  wavelengths.  Attenuation  is  a  function  of  the 
number  of  mean-free-paths  traversed  in  all  cases.  This  approximately 
agrees  with  the  main  features  of  the  UF^  Marteney  emission  data  and 
strongly  suggests  that  UF^  photoabsorption  is  the  common  denominator  be- 
tween the  Florida  and  Marteney  results.  Also  included  in  Appendix  E 
(Figs.  E-5  through  E-8)  are  results  showing  the  original  Marteney  data 
with  the  same  UF^  photoabsorption  strengths  unfolded.  These  unfolded 
data  exhibit  characteristics  in  the  unaltered  arc  emission  data.  However, 
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this  analysis  does  not  resolve  the  disparity  in  emission  coefficient 
magnitudes  as  illustrated  in  Fig.  V-4. 

Shown  in  Fig.  V-4  are  differences  in  the  absolute  values  of  the 
emission  coefficient  measured  at  given  wavelengths.  This  question  has 
been  addressed  by  Schneider,  Campbell,  and  Mack,   where  an  optically 
thin  plasma  and  an  emission  coefficient  that  is  a  direct  function  of 
particle  density  was  assumed.  It  was  then  possible  to  use  the  perfect 
gas  law  to  normalize  the  existing  data  to  a  common  temperature  and  den- 
sity.  (The  exponential  temperature  dependence  of  the  emission  coeffi- 
cient was  neglected.)  The  details  of  this  comparison  are  an  extension 
of  this  study  and  will  only  be  summarized.  The  Marteney  and  Florida 
results  agreed  in  the  1500  to  2200  A  bandwidth.  Beyond  2200  A  the 
Marteney  and  Miller  results  agreed  in  shape  and  magnitude  but  differed 

o 

from  the  Florida  results  by  1  to  3  orders  of  magnitude  at  5000  A. 
However,  if  the  exponential  temperature  dependence  were  considered  in 
the  Miller  (10  000  K)  and  Marteney  (8500  K)  data,  their  magnitudes 
would  differ  as  well.  These  disparities  among  the  three  sets  of  meas- 
urements may,  to  some  extent,  be  attributed  to  UF^  masking  processes. 
Substantial  correlation  in  shape  and  magnitude  resulted  between  Parks' 
theoretical  data  and  the  Florida  experimental  results  when  the  two  were 
normalized  as  previously  described.  The  shape  correlation  can  clearly 
be  seen  in  Fig.  V-4. 

In  summary,  the  emission  coefficient  of  a  uranium  arc  plasma 
(T  '-   8000  K,  P  ■-  0.01  atm)  has  been  measured  and  compared  with  theory 
and  other  similar  experimental  data.  The  Florida  emission  data  were 
measured  from  a  relatively  uncontaminated  uranium  plasma;  the  experi- 
mental comparison  data  were  generated  using  UF^  discharges  exhibiting 
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UF^  photoabsorption  features  that  probably  disguise  the  true  emission 
picture.  The  emission  coefficient  curves  taken  from  the  UF^  discharge 
experiments  have  been,  to  some  extent,  justified  by  relating  their 
shape  to  that  of  the  UF^  photoabsorption  cross  sections.  The  predic- 
tability of  emission  coefficient  shape  and  magnitude  was  addressed  by 
a  comparison  between  the  Florida  results  and  Parks'  calculations. 

o 

Favorable  agreement  exists  in  the  2000  to  4000  A  bandwidth.  Since  no 
calculations  were  made  by  Parks  in  the  vacuum  uv,  this  region  will  be 
investigated  specifically  in  the  next  chapter  by  independent  calcula- 
tions. 


VI.      TUEORETICA L  CALCULA TIONS 
VI- 1 . Introdiicrbion 

Theoretical  justification  for  computing  quantitatively  accurate 
uranium  plasma  emission  coefficients  from  first  principles  is  nearly  im- 
possible because  complete  quantum-mechanical  description  of  the  uranium 
atom  does  not  exist.  There  must  also  be  theory  (statistical  mechanics) 
which  determines  the  level  population  densities  and  ultimately,  emission/ 
absorption  coefficients.  However,  it  is  possible  to  obtain  useful  infor- 
mation about  relatively  strong  emission/absorption  features  as  a  function 
of  wavelength  without  the  approach  described  above.  This  information  is 
in  the  oscillator  strength  calculations  for  those  transitions  in  probable 
and  strong  transition  arrays  located  in  the  wavelength  region(s)  of 
interest.  Oscillator  strength  distribution  can  be  a  reasonable  indicator 
to  the  emission  characteristics  of  the  system  at  those  wavelengths. 

For  singly  ionized  uranium  the  task  of  calculating  transition  array 
oscillator  strengths  is  still  formidable--but  possible.  The  procedure 
used  to  acquire  this  information  for  correlation  to  the  UII  experimental 
emission  coefficients  will  be  discussed. 
VI- 2.      Tei'mi  i  loloni/ 

Consistent  terminology  regarding  atomic  energy  levels  is  necessary 
for  unambiguous  discussions  about  atomic  structure  calculations.  The 
state  of  an  atom  is  the  condition  caused  by  the  collective  motion  of  all 
the  atomic  electrons.  The  state  is  specified  by  four  quantum  numbers  for 
each  electron  or  a  set  of  coupled  quantum  numbers  for  the  entire  atomic 
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system.  The  ground  state  is  the  lowest  energy  state.  A  level  is  repre- 
sented by  the  total  angular  momentum,  J.  The  lowest  energy  level  is 
defined  as  the  ground  level.  Several  states  can  correspond  to  a  given 
energy  level,   rerms  are  collections  of  levels  tagged  by  multiplicity  S 
and  orbital  angular  momentum  L.  The  statistical  weight  (distinct  states) 
in  a  level  is  2J  +  1;  in  a  term  (2S  +  1)(2J  +  1).  Definition  of  the  n  and 
£  quantum  numbers  for  each  electron  orbital  specifies  a  configuration. 
Electrons  in  equivalent  orbital s  are  designated  equivalent  electrons.  A 
transition  of  an  electron  between  two  levels  generates  a  spectral  line, 
whereas  a  multip_let_  is  a  group  of  transitions  between  two  terms.  Finally, 
a  transition  array   is  composed  of  transitions  allowed  between  two  con- 
figurations. 

Coupling  is  the  process  whereby  two  or  more  electron  angular  momenta 
are  combined  into  resultant  angular  momenta.  Regarding  LS,  JJ,  and 
intermediate  coupling,  the  dominance  of  JJ  over  LS  is  expressed  by  the 
relative  magnitudes  of  spin-orbit  and  electrostatic  contributions  to 
energy  separation.  Relative  importance  of  spin-orbit  interaction  gener- 
ally increases  with  increasing  Z  and  n  (principal  quantum  number);  thus, 
for  high-Z  elements  and  large-electron  orbits,  JJ  coupling  would  seem  a 
logical  choice.  Conversely,  for  low-Z  and  small  electron  orbits,  LS 

coupling  would  appear  valid.  However,  at  low-  and  high-Z  there  are  many 

92 
exceptions  to  these  rules  of  thumb;  intermediate  coupling   is  required 

for  many  atomic  systems.  A  good  example  of  the  rule-of-thumb  breakdown 

5 
is  in  the  highest  energy  level  (J  =  2)  in  the  2p  4f  configuration  of 

neutral  neon.  The  electron  eigenfunction  given  in  a  pure  LS  basis  rep- 

93 
resentation  is  composed  of 

0.681|^F2  >  +  0.463|'^D2  >  +  0.567|"'d2  > 
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whereas,  if  pure  LS  coupling  really  existed  (as  expected  of  neon  with 

Z  =  10)  there  would  be  no  contribution  from  two  of  the  three  components 
and  lOO/j  contribution  from  the  third. 
VI- g .    Ccmfi guration  Selection 

The  credibility  of  the  observed  UII  emission  wavelength  can  be  estab- 
lished (to  some  degree)  by  calculating  transition  arrays  and  their  asso- 
ciated oscillator  strength  distribution  of  particular  configuration  pairs. 
Valuable  information  can  be  computed  from  first  principles  relating  to  the 
location,  and  in  some  cases,  the  strength  of  emission  for  a  given  transi- 
tion array.  Generating  specific  transition  arrays  is  directly  related  to 
selection  of  configuration  pairs  that  are  likely  to  produce  spectra  at 
relevant  wavelengths.  The  selection  of  particular  configuration  pairs 
requires  knowledge  of  configuration  average  energies  E  .  The  difference 
AE   between  configurations  is  indicative  of  the  average  transition  array 
wavelength  between  two  specified  configurations.  Average  energies  are 
tabulated  frequently  in  the  literature  for  less  complex  atomic  system  but 
not  for  uranium;  trial  and  error  tactics  were  necessary. 

We  considered  probable  configurations  where  singly  ionized  uranium 

3  2 
could  find  itself,  beginning  with  the  ground-state  configuration  5f  7s 

and  exciting  an  electron  to  another  likely  orbit.  Table  VI-I  summarizes 

configuration  pairs  which  were  ultimately  considered. 

Table  VI-1 
UII  CONFIGURATION  PAIRS 

5f-^7s^  -  5f^7s7p  ^  5f^6d7s^  -  5f^6d^7s 

5f-^7s^  -  5f^7s8p  ^  5f^6d7s^  -  5f^6d7s^ 

5f^7s^  -  5f^6d7s^  5f^6d7s  -  5f-^6d7p  ^ 

5f^7s^  -  5f^7s^7d  5f^7s7p  -  5f'^7s8d  ^ 

5f^6d7s'^  -  5f^6d7s7p  5f^6d7s  -  5f-^7s7p  ^ 


Configuration  pairs  which  were  ^ery   strong  and/or  lead  to  oscillator 
strength  distribution  at  the  desired  wavelengths. 
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Once  defined,  likely  configuration  pairs  were  chosen  (Table  VI-1) 
for  a  quantum  mechanical  calculation  to  establish  the  transition  array 
and  oscillator  strength  distribution  for  a  given  pair.  Appendix  D 
illustrates  the  average  energies  of  these  configurations  as  vjell  as 

average  energies  of  some  lovyer  levels  of  UI  and  UII  as  reported  by 

94 
Radziemski  and  Blaise. 

VI- 4.      Calaul ation  of  Oscillator'  StrenQths 

We  calculated  oscillator  strengths  and  distribution  by  choosing  a 
coupling  representation  to  determine  the  energy  levels  and  wavefunctions 
for  an  atom  in  a  specified  configuration.  For  precision,  energy  levels 
and  wavefunctions  for  all  possible  configurations  must  be  found.  These 
wavefunctions  are  used  to  calculate  the  electric  dipole  matrix  that  leads 
directly  to  oscillator  strength  for  a  given  transition.  The  wavelength 
for  each  transition  is  given  by  the  energy-state  differential  between 
levels  of  the  transition. 

To  determine  energy  states  and  wavefunctions  for  the  atom  (ion)  in 

a  given  configuration  the  Schrodinger  equation,  Eq.  (VI-l)  must  be  solved 

k  k 

for  the  total  wavefunction  ^     and  energy  eigenvalue  E  of  each  state  k. 


H  r 


^^  v^i^ 


(VI-l) 


In  this  formalism  H  represents  the  complete  Hamiltonian  of  the  system. 
This  problem  reduces  to  finding  the  eigenvalues  (E  )  and  eigenvectors 
of  the  Hamiltonian  energy  matrix  in  a  given  basis  representation.  This 
is  accomplished  by  the  diagonal izati on  of  such  a  matrix  into  the  follow- 
ing form: 


<  b|H|b^  >  -  H,,  .  D 
'  '       bb 


(VI-2; 


95 
In  this  case  the  columns  of  the  diagonal izati on  matrix,  D,  form  the 
eigenvectors  and  the  diagonal  elements  are  the  energy  eigenvalues.  The 

standard  "bra"  and  "ket"  notation  used  in  many  quantum  mechanics 

92  95 
texts   '   serves  as  shorthand  for  operator  matrix  elements  H  in  this 

case.  The  proper  definition  of  the  energy  matrix  elements  H,,,  is  the 

bb 

crucial  step  in  obtaining  the  correct  eigenvalues  and  eigenvectors. 
The  complete  Hamiltonian  can  be  written  as 


H  -  -  -"- 
2  m 


Magnetic 

Other 

spin-orbit 

terms 

energy  of 

of 

electron  i 

lesser 

importance 

1=1       i        j>l   ^J    i 

Kinetic       Nuclear  Electrostatic 

energy  of      electro-  potential 

electron       static  energy 

potential  between 

energy  of  electrons 

electron  i  i  and  j 

For  uranium,  Eq.  (VI-1)  cannot  be  solved  in  closed  form  with  this 
Hamiltonian.  Therefore,  approximations  are  usually  made  which  force  a 
given  electron  to  move  in  a  central  field  of  the  nucleus  and  the  N-1 
other  electrons.  If  spherical  symmetry  is  assumed,  the  angular  portions 
of  the  one-electron  wavefunctions  are  hydrogenic  and  can  be  calculated. 
A  self-consistent  field  (SCF)  Hartree-Fock  calculation^^  with  a  spheri- 
cally symmetric  potential  is  used  for  the  radial  wavefunctions.   In  many 
cases  the  Hartree-Fock  approach  is  too  complex  for  solution,  and  further 
approximation  is  necessary.  These  approximations  are  apparent  as  the 
form  of  radial  potential  that  will  account  for  the  system's  exchange 
properties.  It  is  hoped  that  the  potential  used  will  have  the  important 
properties  of  the  Hartree-Fock  potential,  will  yield  "correct"  energy 
eigenvalues,  and  will  establish  properly  orthogonal  determinantal  radial 


III 

i'      wavefunctions 
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When  wavefunctions  and  energy  levels  are  established,  atomic 
spectra  information  can  be  acquired.  Several  types  of  transitions  are 
possible,  but  only  electric  dipole  transitions  will  be  considered  here.. 
Electric  dipole  transitions  are  allowed  only  when  the  matrix  element 
in  the  dipole  matrix  is  nonzero.  The  dipole  matrix  elements  E  are  given 

by 


E  =  <  Vlrlf' 


(VI-4; 


where  >l'  and  T'  represent  the  wavefunction  for  the  upper  and  lower  energy 
state,  and  <  H'|r|H''  >  indicates  the  average  radius  of  an  electron  in  a 
stationary  state  H'  for  the  entire  radiation  process.  The  gf  values 
(statistical  weight  times  oscillator  strength)  are  directly  proportion- 
al to  the  square  of  the  appropriate  dipole  matrix  elements.  The  actual 
coupling  can  affect  significantly  the  strength  and  overall  location  of 
a  given  transition  array  oscillator  strength  distribution.  A  simple 
example  will  illustrate  this  point. 

Consider  the  case  where  the  oscillator  strength  between  the 
■^P.  -^  Sq  levels  is  desired,  where  --^  gf  «  |<  H'lrlH'"  >|  .  Suppose  that 


eigenvector  composition  is 

Level  1  

Level  2  


I  3  I  3 

rP^  >  -  1.0  rP^  > 

I^Sq  >  =  0.9  I^Sq  >  +  0.439  I^Pq  >  . 


*  3  1 

Level  1  is  a  pure   P^  state  and  level  2  is  a  mixture  of  the  S„  and 

3 
Pn  basis  wavefunctions. 


For  a  development  of  the  concept  of  purity,  see  G.  H.  Shortley's 
article,  Ref.  96. 
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Clearly, 


gf  --  I  <  0.9  ^S^  +  0.436  ~^P  Irl  ^P,  >  \'^ 
3     1  0  0  '  '    1    ' 

=  I  0.9  <  ''Sq  |r|  ^P^  >  +  0.436  <  ^P^  \r\    ^P^  >  |^ 

3   I  I  3     ,2 
=  0.19  I  <  Pg  |r|   P^  >  I    . 

In  an  assumed  LS  or  JJ  coupling  scheme,  mixing  would  not  be  mani- 

3     1 
fest  and  the  P.  ->  S^,  vjould  be  disallowed,  resulting  in  gf  =  0. 

Because  the  eigenfunction  for  the  S^^  state  is  a  linear  combination  of 

basis  wavefuiictions,  the  transition  between  levels  1  and  2  is  allowed 

3      3 
by  the  appearance  of  P^  and  P-, .  Mixing  frequently  does  occur  and 

can  play  a  major  role  in  determining  oscillator  strength  distribution 

and  magnitude  of  transition  arrays  between  two  configurations.   For  a 

given  energy  matrix,  transformation  to  JJ,  LS,  or  intermediate  basis 

representation  will  produce  the  same  line  pattern. 

Oscillator  strengths  were  calculated  for  the  UII  configuration 

97 
pairs  labeled  in  Table  VI-1  with  a  computer  code  RCG  created  by  Cowan. 

Briefly,  the  code  implements  the  relativistic  Hartree  exchange  (RHX) 

98 
approach   to  calculate  angular  and  radial  wavefunctions  in  a  basis 

representation  defined  by  the  appropriate  coupling  scheme.   Established 

wavefunctions  and  energy  eigenvalues  are  used  to  compute  wavelengths 

and  magnitudes  of  oscillator  strength  allowed  between  two  configura- 

83 
tions.  Parks  et  al .   attempted  a  similar  calculation  using  a  statis- 
tical model  for  level  distribution  and  the  Thomas- Fermi  approach  for 
the  one-electron  wavefunctions.  Table  VI-2  illustrates  the  major 
considerations  of  the  Parks  relativistic  scaled-Thomas-Fermi  (RSTF) 
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Table  VI-2 
UII  ATOMIC  STRUCTURE  CALCULATIONS  -- 
A  COMPARISON  OF  THE  PARKS  (STF)  AND  RHX  APPROACH 


SUBSTANTIAL  ASPECTS  OF  CALCULATION   PARKS(STF)  RHX 

Scaled-Thonias-Fermi  (STF)  x 

Hartree-Exchange  (HX)  X 

Relativistic  (Dirac)  x       X 

Magnetic  Spin-Orbit  Interaction  X 

Correlation  Effects  X 
Configuration  Interaction 

Self-Consistent  Field  (SCF)  X       X 

Scaled  Thomas-Fermi  Potential  X 

Hartree-Exchange  Potential  x 
Coupling  Schemes 

"I-  ^^one  X  (Statistical 

treatment) 

2.   Intermediate  X 
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0  0 

calculations'   and  the  RHX  calculations.  Results  are  discussed  in 
Sec.  VI-5.  A  flow  diagram  for  RCG  is  given  in  Appendix  B. 

Parks  used  a  relativistic  Hartree-Fock-Slater  self-consistent 
field  treatment  to  calibrate  the  scaled  Thomas-Fermi  potential. ^^  With 
this  potential,  he  computes  one-electron  wavefunctions  for  the  opacity 
calculations.  The  Thomas-Fermi  model  for  heavy  elements  is  adequate 
until  encountering  electrons  in  overlapping  levels  and  when  matching 
large  and  small  radii  boundary  conditions.  The  scaled  Thomas-Fermi 
potential  does  not  account  for  the  exchange  contribution  to  the  electron 
potential  energy.  The  Hartree-exchange  potential  resembles  the  charac- 
teristics of  the  Hartree-Fock  potential,  and  accounts  for  the  exchange 
terms  in  the  central-field  approximation.  This  implies  that  wavefunc- 
tions developed  with  the  Thomas-Fermi  potential  often  will  exhibit  a 
greater  uncertainty  than  those  calculated  with  the  Hartree-exchange  po- 
tential. Even  the  energy  eigenvalue  calculations  can  be  poor,  as  shown 
for  the  5f  uranium  binding  energy  in  Fig.  1  of  Ref.  83.  The  exclusion 
of  correlation  effects  in  the  Parks  model  is  not  a  serious  drawback  and 

no 

may  be  inconsequential;   however,  omission  of  a  detailed  treatment  of 
actual  coupling  could  be  its  most  serious  shortcoming.  Large  inaccura- 
cies are  possible  in  line-pattern  calculations  and  oscillator  strengths, 
particularly  if  one  configuration  in  the  transition  array  is  not  well 
characterized  by  a  pure  coupling  scheme.  This  situation  arises  often 
in  the  case  of  uranium. 
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VI- 3.      Compan-son  of  Results 

Figure  VI-1  is  a  collective  plot  of  the  RHX  oscillator  strength 
^  distributions,  the  Parks  theoretical  emission  coefficient  STF  calcula- 

0  0 

tion,   and  the  experimentally  determined  UII  emission  coefficient,  all 
as  a  function  of  wavelength.  Only  the  RHX  oscillator  strengths  relate 
to  the  absolute  values  found  on  the  ordinate,  whereas  the  Florida  and 
Parks  emission  coefficient  curves  are  slightly  distorted  in  the  rela- 
tive magnitude  but  the  major  inflection  points  have  been  preserved. 
This  plot  format  shows  only  the  significant  features  of  the  experimen- 
tal emission  coefficient  curves  for  comparison  to  theory;  magnitude 
was  discussed  in  the  previous  chapter.  For  each  transition  array 
oscillator  strength  distribution  (TAOSD)  shown  in  Fig.  VI-1,  the  TAOSD 
and  the  corresponding  fractional  number  of  spectral  lines  per  wavelength 

interval  are  plotted  in  Appendix  C.  Plots  showing  the  bf^ls^  -   Bf-^/sSp 

3  2     3 
and  5f  7s  -  5f  7s7p  transition  arrays  clearly  show  that  many  spectral 

lines  in  a  cluster  do  not  insure  strong  oscillator  strength  and  strong 
emission. 

In  general,  the  Parks  and  RHX  calculations  support  the  variation  in 
shape  in  the  Florida  experimental  emission  coefficient.  Several  posi- 
tive points  are  made  when  comparing  the  RHX  calculations  with  the 
Florida  measurements. 

(1)  For  our  configurations,  the  strongest  emission  is  likely  to 
occur  in  the  2500  to  5000  A  bandwidth. 


(2)  An  overall  peak  is  likely  in  TAOSD  at  2700  "^  ^|?J? 


relatively  smaller  peak  at  2041  "^  ^11   A. 
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(3)  The  expey^iiiiental  emission  coefficient  "shelf"  between  1900  and 

°  3 

1750  A  can  bo  supported  by  the  fairly  strong  TAOSD  from  the  5f  8s7p  - 

3 
5f  7s8d  configuration  pair. 

o 

(4)  The  emission  coefficient  decrease  at  ''-^   2500  A  is  validated  by 

o 

a  corresponding  sparcity  of  oscillator  strength  between  2200  and  2500  A. 

(5)  Similarly,  an  observed  emission  coefficient  decrease  in  the 

o 

1600  to  1700  A  bandwidth  is  supported  by  a  lack  of  oscillator  strength 
with  similar  qualification  as  mentioned  in  point  4. 

(6)  The  5f^7s^  -  5f'^7s8p  TAOSD  indicates  good  potential  for  a  small 

+  74  ° 
but  significant  emission  peak  at  -  1539  _  go  A,  and  it  represents  a 

theoretical-experimetnal  verification  of  prominent  UII  emission  in  the 
vacuum  ultraviolet  wavelength  region. 

The  Parks  calculation,  as  far  as  it  was  taken  in  wavelength,  re- 
produces the  experimental  results  at  the  major  peaks.  He  attributes 

o 

the  appearance  of  distinct  emission/absorption  peaks  at  2800  and  2200  A 

3  2     3 

to  the  transition  array  formed  between  the  5f  7s  -  5f  7s7p  configuration- 

83 
pair   for  J  =  1/2,  3/2  of  the  7p  electron.  One  can  hypothesize  a 

corresponding  enhancement  of  the  oscillator  strength  at  these  wavelengths 

for  this  transition  array.  However,  the  RHX  calculation  for  the  same 

transition  array  shows  three  clusters  of  lines  (Appendix  C,  Fig.  C-10) 

at  5890  A,  4000  A,  and  2700  A,  but  only  one  distinct  oscillator  strength 

O 

peak  at  2700  A  (Appendix  C,  Fig.  C-9).   It  is  evident  that  the  distinct 
peaks  predicted  by  the  Parks  model  for  J  =  1/2,  3/2  of  the  7p  electron 
are  not  predicted  by  the  more  sophisticated  model  of  Cowan.  The  experi- 

o 

mental    peak  at  2100  A  is   predicted  by  the  RHX  calculation  of  the  transi- 

3  3 

tion  array  formed   from  the  5f  7s7p  -   5f  7s8d  configuration-pai r. 
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The  disparity  in  peak  identification  between  the  Parks  and  Cowan 
models  lies  in  the  theoretical  approach.   Emission  peak  location  is 
determined  by  level  splittings  caused  principally  by  spin-orbit  and 
electrostatic  (direct  and  exchange)  interactions.   For  UII  Parks  ne- 
glects the  exchange  portion  of  the  electrostatic  interaction  energy 
altogether;  he  also  neglects  the  higher-order  direct  contributions  as- 
sociated with  equivalent  electrons  in  the  5f  orbital.  The  RHX  calcu- 
lation indicates  that  not  only  are  these  interaction  energies  impor- 
tant, they  are  dominant,  particularly  with  regard  to  transition  array 
line-cluster  location  and  oscillator  strength  distribution.  The  RHX 
coupled  wavefunctions  are  certainly  more  accurate  than  those  uncoupled 
wavefunctions  used  in  the  Parks  treatment.  Most  quantities  dependent 
upon  wavefunction  description  (including  oscillator  strengths)  reflect 
the  accuracy  (or  inaccuracy)  of  the  wavefunctions  used.   Further,  Parks 
includes  level  populations  in  his  calculations,  but  the  main  effect 
of  ignoring  this  feature  is  to  produce  changes  in  relative  peak  magni- 
tude when  calculating  emission  coefficients,  but  not  in  the  locations 
of  such  peaks.  However,  the  Parks  calculations  indicate  strong  emission 

o 

trends  between  2300  and  4000  A,  and  in  that  sense,  are  supported  by  ex- 
periment and  RHX  theory.   In  fact,  it  is  remarkable  that,  for  the  most 
part,  the  Parks  calculations  came  out  as  well  as  they  did  considering 
the  approximations  used  and  the  complexity  of  the  problem. 


VTI.    CLOSING  REMARKS 

Emission  coefficients  of  uranium  plasma  at  two  distinct  tempera- 
tures and  pressures  over  an  unusually  broad  wavelength  region  have  been 

measured.  These  new  data  represent  an  extension  of  previous  information 

79       22  21       Rl 

published  by  Krascella,   Miller,   Marteney  et  al . ,   Randol  ,   and 

CO 

Kylstra.    The  results  are  significant  additions  to  available  informa- 
tion, and  the  major  points  are  summarized  below. 

Specific  emission  coefficients  for  a  low-pressure  uranium  arc  at 

o 

5500  K,  0.001-atm  uranium  pressure  were  measured  from  2500  to  5500  A. 

These  data  were  in  reasonable  agreement  with  the  semi-empirical  work 

79 
of  Krascella   and  the  spectral  line  catalogue  of  the  first  spectrum 

of  uranium  reported  by  Steinhaus  et  al.   '    For  this  arc  there  was  no 

o 

emission  detected  below  2500  A,  probably  because  of  the  relatively  low 
temperature  and  the  importance  of  radiation  from  neutral  uranium.  Plasma 
temperature,  pressure,  and  emission  coefficient  were  obtained  simultan- 
eously, thus  reducing  fluctuation  problems  usually  associated  with 
photoelectric  diagnostics.  The  temperature  and  pressure  diagnostics 
depended  upon  assuming  LTE,  and  although  not  proven,  LTE  was  reasonably 
substantiated  with  cross-checks  and  comparison  to  NBS  work  by  Voigt. 

Specific  etnission  coefficients  of  a  high-pressure  helium-uranium 
arc  at  approximately  8000  K,  0.01-atm  uranium  pressure  were  measured 

o 

from  1050  to  6000  A.  These  results  represent  a  rare  successful  attempt 
to  obtain  emission  coefficient  data  for  singly  ionized  uranium  in  the 
vacuum  ultraviolet.   Intensity  calibration  in  the  visible  and  vacuum 

104 


105 


ultraviolet  v«s  done  with  a  tungsten  standard  and  hydrogen  discharge  for 
each  respective  wavelength  region.  Calibration  below  the  wavelength 
cutoff  (1750  A)  of  the  hydrogen  standard  was  inferred  by  extrapolation. 
Direct  comparisons  with  similar  experimental  results  by  Marteney  et  al."^ 
were  understandably  poor  because  of  differing  discharge  systems  and  plasma 
characteristics.  The  Marteney  data  were  consistent  in  shape  with  our 
results  because  of  probable  UFg  photoabsorption  in  the  Marteney  experi- 
ment. The  emission  coefficient  values  of  the  Florida  and  Marteney  efforts 
agreed  to  within  an  order  of  magnitude  if  scaled  to  a  common  temperature 

A    ^  -f   91 

and  densi ty. 

Comparison  of  our  experimental  results  for  singly  ionized  uranium 

0  0 

with  Parks  et  al."  theoretical  calculations  using  a  relativistic  scaled 
Thomas-Fermi  model,  and  those  using  the  relativistic  Hartree  exchange 
approach,  provided  insight  into  the  origin  of  the  vacuum  uv  emission  and 
strong  emission  at  other  wavelengths.  The  Parks  calculation  included  a 
statistical  mechanical  treatment  of  the  calculated  energy  level  popula- 
tion. The  RHX  computations  defined  only  the  oscillator  strength  distri- 
bution for  selected  transition  arrays,  which  are  useful  for  predicting 
significant  emission  location  as  a  function  of  wavelength.  The  Parks 
and  RHX  calculations  supported  the  experimental  results  to  varying  degrees 

o 

The  peak  locations  at  2100  and  2900  A  were  predicted  by  both  approaches, 
but  for  different  reasons.  The  disagreement  is  explained  by  the 
distinctions  between  the  two  models.  Although  the  Parks  calculations 
did  not  extend  through  the  vacuum  uv,  RHX  calculations  conclusively 

o 

predict  the  peak  in  UII  emission  coefficient  observed  at  about  1540  A. 
These  theoretical-experimental  comparisons  indicate  that  substantial 
success  can  be  anticipated  in  predictions  of  qualitative  features  in  the 
wavelength  dependence  of  emission  from  plasmas  of  very  complex  systems. 
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This  effort  has  resulted  in  a  credible  definition  of  uranium  plasma 
emission  properties  over  an  extended  wavelength  range  near  the  reported 
temperature/pressure  regimes.  Presumably,  the  success  of  the  theoretical 
predictions  will  increase  the  confidence  in  atomic  structure  calculations 
for  heavy  elements  and  will  indicate  possible  areas  for  model  improvement. 
Although  there  were  numerous  difficulties  in  this  research,  and  many 
compromises  were  made,  the  overall  result  is  clearly  progressive  in  the 
experimental  and  theoretical  aspects  of  the  problem. 


APPENDIX  A 
SAHA   NUMBER  DENSITIES  AND  NORMAL  TEMPERATURES 

Saha  analysis  was  used  to  calculate  LTE  number  densities  of 
uranium  for  two  stages  of  ionization  at  several  total  pressures  of 
interest.  The  Saha  equations  are  given  by 

S.(T)---^- g.-^ T    e  T     ,      (A-1) 

1       1        h 


where 


i  =  the  ionization  stage, 

-3 


n^   =  number  density   (cm     )   for  all    i-fold  ionized  particles, 


3, 


n  -   electron  number  density  (cm~  ) 


U^.  (T)  -  partition  function  of  i-fold  ionized  atom, 

-  ?s 
m^  ==  electron  rest  mass  (=  9.108  x  10    gm), 

h  ^  Planck's  constant  (=  6.626  x  10    ergs-seconds), 

k  -  Boltzmann's  Constant  (=  1.381  x  10"^^  ergs/K), 

E^  =  ionization  energy  for  ionization  from  i  ->-  i+1, 

T  =  absolute  temperature  (K). 

The  number  density  solutions  were  obtained  using  Eq.  (A-1)  coupled  with 

the  equation  of  charge  neutrality: 

i 

^e  =  ^  """i   '  (A-2) 

i-0   ^ 

and  the  equation  of  state: 

i 

P  =  kT  {n  +  E  n.}  .  (A-3) 

i  =  0  ^ 
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These  number  densities  are  plotted  in  Figs.  A-1  through  A-6  as  a  func- 
tion of  plasma  temperature.   Uranium  ionization  potentials  used  were 
those  reported  by  Williamson  et  al.    and  partition  function  tempera- 
ture dependence  was  included  by  incorporating  semiempi rical  curves 
defined  by  Krascella  et  al.    Lowering  the  ionization  potentials  was 
not  considered. 

The  Saha  number  densities  at  various  total  pressures  were  then 
used  to  define  normal  temperatures  for  representative  UI  and  UII  transi- 
tions. The  mathematical  formalism  is  discussed  in  Sees.  II-4B  and  4C-- 
this  information  is  in  Fig.  A-7. 
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Fig.    A-1.       Uranium    plasma    Saha    number    densities 
P  =  3    a  till. 
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Fig.    A -2.       Uranium    plasma    Saha    number    densities 
P=0.3    atm. 
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Fig.  A-3.   Uranium  plasma  Saha  number  densities, 
P=0.1  atm. 
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Fig.  A-4.   Uranium  plasma  Saha  number  densities 
P  =  0.01  atm. 
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Fiq.  A-5.   Uranium  plasma  Saha  number  densities, 
P=0.001  atm. 
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Fig.  A-6.   Uranium  plasma  Saha  number  densities 
P-0.0001  atm. 
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Fig.  A-7.   Pressure  vs  normal  temperature  of 

neutral  and  singly- ionized  uranium  in  LTE 


APPENDTX  B 
COMPUTER  SCHEMATIC  FOR  THE  COWAN  RCG   CODE 
The  information  in  Fig.  B-1  was  taken  from  two  sources^^'^°°  and 
represents  the  major  steps  in  the  RCG  code  for  calculating  oscillator 
strength  distributions  used  in  this  study.   It  is  not  meant  to  expose 
the  intricacies  of  the  entire  calculation  but  to  present  an  overview  of 
the  calculations  reported  in  this  study.  Electrostatic  interaction 
parameters  in  the  form  of  Slater  radial  integrals^^  are  denoted  by  f'^ 
and  G  for  the  direct  and  exchange  contributions,  respectively.  The 
spin-orbit  interaction  term  is  indicated  by  '£,   and  E^   is  the  average 
energy  of  all  states  of  a  configuration. 
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Code  schematic  for  RCG  calculation  of 
oscillator  strength  distributions. 


APPENDIX  C 
SPECTPAL  LINE  AND  OSCILLATOR  STRENGTH  DPSTRIBUTIONS 
I'VE  SELECTED  UII  CONFIGURATION  PAIRS 

The  following  figures  illustrate  UII  spectral  line  and  oscillator 
strength  distributions  as  a  function  of  wavenumber.  They  compose  the 
majority  of  graphical  output  from  the  RCG  calculation.  Figures  C-1  , 
C-3,  C-5,  C-7,  and  C-9  are  the  oscillator  strength  distributions  for  a 
specified  UII  configuration  pair  and  have  already  been  presented  in 
Fig.  VI-2.  Figures  C-2,  C-4,  C-6,  C-8,  and  C-10  illustrate  the  calculated 
fractional  spectral  line  distributions  for  corresponding  configuration 
pairs.  The  number  of  spectral  lines  calculated  for  the  five  configura- 
tion pairs  reported  approached  203  000.  Only  those  lines  which  had 
strengths  above  a  user-set  cutoff  were  plotted. 
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Fig.  C-1.   RHX  transition  array  oscillator  strength 
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APPENDIX  D 

HEPRESENTATTVE  ENERGIES  OF  SEIECTED 

VI  AND  UII  CONFIGURATIONS 

When  estimating  probable  wavelength  locations  of  transitions 

,;| 

I       between  configut^ation  pairs,  the  difference  in  average  configuration 

li! 

energies  is  useful.  Radzieniski  and  Blaise  ^   have  tabulated  energies  of 
the  lowest  known  levels  of  ten  configurations  of  UI  and  UII  and  this 
information  is  shown  in  Figs.  D-1  and  D-2.  While  the  lowest  level  energy 
is  not  the  true  average  configuration  energy,  it  is  usually  a  good  approx- 
imation. Also  in  Fig.  D-3  are  the  average  configuration  energies,  cal- 
culated by  the  RHX  method,  of  the  UII  configurations  used  in  this  study. 
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APPENDIX  E 

UF,.  PHOTOABSOHPTION  EFFECTS  OF  THE  MARTENEY  AND 
0 

FLORIDA   EMISSION  COEFFICIENT  DATA 


UFg  photoabsorption  strengths  as  a  function  of  wavelength  were 
calculated  using  the  cross-section  data  of  DePoorter,   ' 
Srivastava,  ^  McDiarmid,^''  Trajmar,^^  and  Hay^^  (Figs.  V-5  and  V-6). 
The  intensity  of  transmitted  light  I-^  that  has  passed  through  a  de- 
fined depth  of  UFg  of  given  molecular  concentration  is  related  to  the 
incident  light  Iq  by  the  equation 


where        £  =  UFg  cell  or  layer  thickness  in  cm, 

3 

n  =  molecular  concentration  in  UF^  molecules/cm  ,  and 

6 

o{\)   =  v/avelenght  dependent  UFg  photoabsorption  cross- 

2 
section   in  cm  . 

These  absorption  factors  (Ij/Iq)  were  used  to  illustrate  the  possible 

effect  of  UFg  photoabsorption  on  the  arc  emission  coefficient  data  of  the 

present  experiment.  This  was  done  by  folding  absorption  strenghts 

calculated  using  Eq.  (E-1)  with  the  arc  emission  coefficient  values  at 

regularly  defined  wavelengths.  The  folding  was  preformed  at  UF^  layer 

thicknesses  of  1.5,  2.0,  3.0,  and  5.0  cm  with  molecular  concentrations 

ranging  from  1  x  10   -^  1  x  10   cm"  .  The  results  of  these  calculations 

are  shown  in  Figs.  E-1  through  E-8. 
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134 
Presumably,  the  Marteney  data  have  UFg  photoabsorption  contamina- 
tion at  the  strenghs  defined  by  UFg  layers  of  1.5  to  5.0  cm  and  molecular 
concentrations  of  1  x  10^^  to  1  x  10^''  cm'^.  The  reciprocal  absorption 
strength  factors  were  also  folded  into  the  matreney  data  to  remove  the 
potential  UFg  photoabsorption  effects,  thus  yielding  an  uncontaminated 
uranium  emission  coeffienent.  These  results  are  displayed  in  Figs.  E-5 
through  E-8. 
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